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BROADBAND  MODELS  FOR  PREDICTING  BISTATIC  BOTTOM, 
SURFACE,  AND  VOLUME  SCATTERING  STRENGTHS* * 


1.  INTRODUCTION 

For  a  low-frequency  (50  to  1000  Hz)  or  mid- frequency  (1  to  10  kHz)  operational  active  sonar,  scattering 
from  the  ocean  boundaries  and  biologies,  coupled  with  propagation  conditions,  can  severely  limit  the  detect¬ 
ability  of  returns  from  features  of  interest.  Furthermore,  reverberation  levels  can  vary  dramatically,  depending 
on  the  local  geology,  oceanography,  and  biology  Hence,  making  accurate  predictions  of  active  sonar  perfor¬ 
mance  will,  in  turn,  depend  on  finding  suitable  empirical  or  physics-based  models  that  accurately  describe  the 
scattering.  This  paper  presents  broadband,  three-dimensional  (3D),  ocean  bottom,  surface,  and  volume  (fish) 
scattering  models  that  should  not  only  improve  the  robustness  and  accuracy  of  sonar  performance  predictions  of 
multistatic  reverberation,  but  also  provide  new  insights  as  to  the  significance  of  each  scattering  mechanism. 

Bottom  scatter  is  generally  the  dominant  interference  source  for  active  systems  in  littoral  waters.  However, 
it  is  a  poorly  understood  combination  of  scattering  from  rough  interfaces  (water-sediment,  sediment-sediment, 
sediment-basement)  and  the  sediment  volume.  This  is  a  result  of  both  the  complexity  of  the  ocean  bottom  and 
the  difficulty  of  acquiring  in  situ  environmental  inputs  to  the  accuracy  required  by  models. 

Surface  scattering  is  caused  by  the  interaction  of  acoustic  energy  with  environmental  features  at  or  near  the 
ocean  surface.  The  dynamic  nature  of  the  air-sea  boundary  interaction  zone  complicates  this  process  [1].  When 
wave  breaking  is  significant,  air  becomes  entrained  in  the  form  of  subsurface  bubbles.  Under  these  conditions, 
both  the  rough  air-sea  interface  and  bubble  clouds  may  contribute  to  the  acoustic  scattering. 

Volume  scattering  at  low-  and  mid-frequencies  is  primarily  caused  by  acoustic  interaction  with  fish 
swimbladders.  General  scattering  models  for  dispersed  bladdered  fish  (incoherent  sum  of  scattering  from 
individual  fish)  have  been  developed  [2]  that,  when  coupled  with  fish-distribution  algorithms,  provide  regional 
predictions  of  mean  volume  scattering  [3].  (Fish  distribution  and  behavior  are  linked  to  oceanographic  variables 
such  as  ocean  temperature,  mixed-layer  depths  and  topography  that  can  be  used  to  advantage  to  predict  fish 
distributions.)  When  fish  school,  as  is  common  in  shallow  water,  the  scattering  picture  is  more  complicated. 
The  Naval  Research  Laboratory  (NRL)  has  developed  a  scattering  model  for  small-sized  fish  schools  that 
includes  all  orders  of  multiple  scattering  interactions  between  the  fish  (with  the  aggregate  scattering  field 
calculated  by  coherent  summation)  [4] .  Another  set  of  complications  arises  when  fish  are  in  the  presence  of 
reflecting  boundaries  [5].  This  paper  extends  the  free-field  dispersed- fish  model  by  incorporating  boundary- 
interference  effects  for  both  the  ocean  surface  and  bottom. 

Recent  broadband  scattering  measurements  conducted  by  NRL  under  the  Space  and  Naval  Warfare  Systems 
Command’s  Critical  Sea  Test  (CST)  program  [6]  (1988-1996)  and  the  Office  of  Naval  Research’s  Littoral 
Warfare  Advanced  Development  (LWAD)  project  [7]  (1996-present),  have  helped  to  distinguish  the  important 
scattering  mechanisms  and  determine  the  environmental  conditions  for  which  each  mechanism  makes  signifi¬ 
cant  contributions.  Using  a  variety  of  systems  and  waveforms,  NRL  made  direct-path  measurements  of  bottom, 
surface,  and  volume  backscattering  characteristics  for  a  range  of  frequencies  and  environmental  conditions  in 
both  the  open  ocean  and  shallow  water.  These  measured  acoustic  quantities  have  been  correlated  to  archival 
bottom  properties,  environmental  measurements  of  the  air-sea  boundary  conditions  (such  as  wind  speed,  wave 
spectra,  and  bubble  parameters),  and  in  situ  assessments  of  biological  quantities  (such  as  fish  densities  and 
sizes). 
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In  conjunction  with  these  scattering  measurements,  NRL  has  been  developing  both  empirical  and  physics- 
based  models  of  scattering  strength,  with  an  emphasis  on  generating  formulas  suitable  for  inclusion  as  submodels 
in  both  active  performance/reverberation  models  and  inverse  algorithms  (such  as  environmentally  adaptive 
techniques).  To  this  end,  this  paper  describes  the  latest  versions  of  these  models  for  estimating  the  dependence 
of  scattering  strength  on  the  incident  and  scattered  grazing  angles,  the  bistatic  angle,  the  acoustic  frequency, 
biological  descriptors  of  the  fish,  and  physical  descriptors  of  the  environment.  After  a  few  definitions,  the  paper 
begins  with  a  discussion  of  3D  interface-scattering  models,  applicable  to  both  ocean  boundaries.  Then,  dis¬ 
cussed  in  turn,  are  the  new  models  for  calculating  the  scattering  strengths  of  the  bottom  interface;  dispersed 
bladdered  fish  near  boundaries;  and  the  sea  surface  (air-sea  interface  and  bubble  clouds).  We  conclude  with  a 
few  comments  and  recommendations. 

2.  DEFINITIONS 

Scattering  strength  (in  dB)  is  defined  by 


£S  =  HMog10a,  (1) 

where  a  is  the  scattering  cross  section  per  unit  area  (per  unit  solid  angle),  a  dimensionless  quantity.  Scatter¬ 
ing  strength  represents  a  mean  far-field  (plane- wave)  quantity. 

The  scattering  geometry  adopted  in  this  paper  follows  Jackson  [8]  and  is  illustrated  in  Fig.  1  for  the  ocean 
bottom.  (For  the  ocean  surface,  the  geometry  is  inverted,  so  that  positive  z  corresponds  to  depth  below  the 
surface.)  Here,  k  and  q  are  the  incoming  and  outgoing  wavevectors,  respectively,  and  define  the  incident  grazing 
angle  Qinc  and  the  scattered  grazing  angle  Qscat.  The  bistatic  angle  §bi  is  defined  as  the  difference  in  azimuth 
between  the  incident  and  scattered  directions.  (There  is  only  one  bistatic  angle;  the  scattering  statistics  of  the 
interface  roughness  and  volume  heterogeneities  are  assumed  to  be  locally  transversely  isotropic.)  This  in  turn 
implies  SS  will  be  an  even  function  of  §bi,  i.e.,  SS(Qinc,  Qscat,  §bi  )  =  SS(Qinc,  Qscat,  -(),„)• 

With  these  definitions,  we  have  0  <  Qinc  <  90°,  0  <  Qscat  <  180°,  and  -180°  <  §bi  <  180°.  Bistatic  angles 
of  ±180°  correspond  to  the  backscatter  direction,  and  a  bistatic  angle  of  0°  with  Qinc  =  Qscat  corresponds  to  the 
specular  direction.  (Figure  1  also  shows  that  SS(Qinc,  Qscat,  <\>bi)  =  SS( Qinc,  180°  -  Qscat,  <| >w  ±  180°),  so  that  <\>hi 
=  ±180°  with  Qscat  =  180°  -  Qinc  also  correspond  to  specular.)  All  the  scattering  strength  models  will  satisfy 
reciprocity  in  that  SS(Qinc,  Qscat,  <$>bi)  =  SS(Qscat,  Qinc, 

In  terms  of  angles,  the  incident  and  scattered  wavevectors  can  be  written  as 

k  =  k0(cosdincx-sinQincz),  (2) 

q  =  k0  (cos0JC<tf  cos  tybix  +  cos  Qscat  sin  tybiy  +  sin  dscatz),  (3) 

where  k0  =  2nflc0  is  the  acoustic  wavenumber,  with  /  the  acoustic  frequency  and  c0  the  speed  of  sound  in 
bubble-free  seawater  at  the  interface.  Later  expressions  for  the  interface  scattering  cross  sections  will  depend  on 
the  horizontal  and  vertical  magnitudes  of  the  difference  of  these  two  wavevectors.  Defining  k  -  q  =  (Qh,  Q:), 
their  magnitudes  can  be  expressed  in  terms  of  the  scattering  angles: 


\Qh\  =  k0^[cos^ 


®inc  +  COS 


-2cos0,„ccos0Jca;  COS  <\>bi, 


(4) 


Q-  =-^o(sin0(„c+sin0SCflf). 


(5) 


MKS  units  are  used  throughout  this  paper. 
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SIDE  VIEW 


Fig.  1  —  Definition  of  scattering  angles:  3D,  side, 
and  top  views,  as  illustrated  for  the  ocean  bottom 


3.  INTERFACE  SCATTERING  MODELS 

A  variety  of  models  have  been  applied  to  sound  scattering  by  rough  ocean  boundaries,  such  as  first-order 
(Born  approximation)  and  higher-order  perturbation  theories,  and  composite-roughness  models.  For  this  paper, 
the  small  slope  formalism  of  Voronovich  [9-11],  as  adapted  for  the  ocean  surface  by  Dashen  et  al.  [12]  and 
elastic  seafloors  by  Wurmser  [13],  is  used  to  predict  the  contributions  of  the  boundary  interfaces  to  the  total 
scattering.  This  lowest-order  small-slope  approximation  (SSA)  models  interface  scattering  strength  to  all  orders 
in  surface  height  h,  and  through  first  order  in  the  derivatives  of  h  (surface  slope).  Using  the  local  small-slope 
approximation  instead  of  the  standard  first-order  perturbation  approximation  increases  the  accuracy  of  the 
prediction  at  the  expense  of  a  modest  increase  in  numerical  complexity  [14]. 

Bottom  and  surface  roughness  are  taken  to  be  Gaussian  random  processes,  each  described  by  a  two- 
parameter,  isotropic,  2D  roughness  spectral  density  of  the  form  [8,12] 


W(Qh)  = 


^2 

(*b| Qh\y2  ’ 


(6) 


with  y2  e(2,  4),  \Qh\  given  by  Eq.  (4),  and  h0  a  normalizing  reference  distance  of  1  m.  (Each  subscript  “2” 
refers  to  the  boundary  relief  being  described  by  a  2D  random  process  [8].)  The  units  of  w2  are  m4.  Higher 
values  of  y2  in  the  given  range  correspond  to  a  weaker  frequency  dependence. 
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As  is  standard,  the  coherent  specular  component  is  removed  in  this  treatment  of  rough  surface  scattering 
theory.  Following  Jackson  [8],  these  interface-scattering  models  assume  that  all  scattered  energy  is  incoherent. 
Hence,  it  would  be  a  mistake  to  add  a  coherent  reflected  wave  to  these  models,  as  conservation  of  energy  would 
then  be  violated.  This  is  an  artifact  of  the  spectrum  given  in  Eq.  (6).  For  such  spectra,  there  is  no  small- 
wavenumber  (i.e.,  large-scale)  cut-off,  and  so  the  surface  remains  rough  even  at  the  largest  scales.  Specular 
reflection  is  ordinarily  generated  because  the  surface  flattens  out  at  the  largest  scales.  On  the  other  hand,  if  the 
surface  has  very  little  large-scale  roughness,  then  there  will  be  a  sharp  peak  in  aint  in  the  specular  direction, 
with  its  magnitude  and  width  such  as  to  give  the  correct  energy  flux.  Hence,  issues  related  to  the  large  scales 
(i.e.,  small  wavenumbers),  most  notably  the  lack  of  a  low- wavenumber  cut-off,  generally  affect  scattering  near 
the  specular  direction.  (We  are  working  on  improvements  to  these  models,  but  a  fundamental  problem  remains 
with  the  ocean  bottom  since  there  is  no  general  procedure  for  determining  the  cut-off.) 

Because  the  model  development  is  applicable  to  both  the  ocean  bottom  and  surface,  we  begin  with  the 
general,  lowest-order  SSA  formulation  and  discuss  its  numerical  implementation  before  specializing  to  specific 
boundaries.  These  developments  are  applicable  to  the  full  3D  problem. 

Theoretical  Formulation 


The  SSA  result  for  the  noncoherent  component  of  the  scattering  cross  section  per  unit  area  (per  unit  solid 
angle)  for  a  random,  rough  interface  has  the  form  [12] 


1 

P 

8rc 

\Qh\Qz 

(7) 


where  (3  is  an  algebraic  form  that  depends  on  the  boundary  conditions  prevailing  at  the  interface,  I  is  an  integral 
involving  the  spatial  spectrum  of  the  roughness,  and  \Qh\  and  Qz  are  given  by  Eqs.  (4)  and  (5),  respectively. 
(Equation  (7)  is  similar  in  form  to  Eq.  (38)  of  Jackson  et  al.  [15]  and  Eq.  (11)  of  Williams  and  Jackson  [16], 
except  that  the  coefficient  in  front  of  the  integral  is  a  small-slope  approximation  while  theirs  are  Kirchhoff 
approximations.) 

With  v  =  (y2  -  2)/2,  and  thus  v  e(0,  1),  the  integral  in  Eq.  (7)  is  [17-18] 


1(a)  =  JQ  J0 (y)yB(y, a )dy,  where  B(y, a)  =  exp(-ay2v) , 


(8) 


/0  is  the  zeroth-order  Bessel  function  of  the  first  kind,  and  a  is  given  by 

\2 


a  = 


{KnsQz)  rq-v) 

(2/*0|e„|)2v  r(i+v>  ’ 


with  the  mean-square  roughness  /i,2,ilv  defined  by 


2  1 W2_ 

nrms  1 2 


hz0v 


(9) 


(10) 


(In  the  derivation  of  the  lowest-order  SSA,  a  finite  cut-off  spectrum  is  used.  However,  as  a  result  of  cancella¬ 
tions  during  the  procedure,  the  finite  cut-off  disappears  and  W  effectively  assumes  the  pure  power  law  form 
described  by  Eq.  (6).) 

Assuming  the  sound  speeds— both  the  real  and  imaginary  (attenuation)  parts— are  frequency-independent, 
then  the  only  frequency  dependence  in  aint  comes  from  the  parameter  a,  which  under  these  circumstances  is 
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proportional  to  /2(1  v)  =  /4  72  [18].  In  this  case,  the  only  environmental  parameter  that  influences  the  fre¬ 
quency  dependence  is  the  boundary  roughness  spectral  exponent  y2.  On  the  other  hand,  once  the  real  or 
imaginary  part  of  a  sound  speed  acquires  a  dependence  on  frequency,  then  a,  (3,  and  of  course  aint  all  acquire 
complicated  frequency  dependencies. 

Numerical  Implementation 


While  the  integral  in  Eq.  (7)  is  well-defined,  it  is  not  easily  evaluated  due  to  an  oscillating  integrand  that 
decays  at  a  highly  variable  rate.  This  difficulty  is  compounded  by  the  requirement  that  1(a)  be  evaluated  over  a 
wide  range  of  a  (a  wide  range  of  grazing  angles).  Two  approaches  have  been  developed  to  numerically 
implement  Eq.  (7).  Together,  they  provide  a  recipe  for  calculating  aint  that  covers  most  y2  of  practical  interest. 

Two-Series  Approach 


The  first  approach  [19]  uses  a  series  representation  of  J0(y)  to  perform  the  integration  of  1(a)  term-wise  for 
large  a,  and  a  second  series  for  B(y ,  a)  to  do  the  same  for  small  a.  (The  power  law  form  of  the  argument  of  the 
exponential  in  B(y ,  a)  occurs  because  the  autocorrelation  function  for  the  surface  height  function  has  been 
expanded  to  first  order.  Numerically,  this  form  is  particularly  useful  because  at  this  order,  the  finite  cut-off  of 
the  surface  spectrum  disappears,  and  because  there  are  no  general  rules  for  specifying  the  spectral  peak  for  the 
ocean  bottom.)  With  this  approach,  the  integral  of  Eq.  (8)  is  approximated  as  a  pair  of  easy-to-evaluate  sums 
[17-19]. 


1(a) 


■  2yr  y.+i  r(»/v) 
v„=i  r2(n) 


(4a1/v)" 


a  <a 


(Ha) 


~  —  (~)”+1  sin(7mv)  —  +  —  (4va)  a<a.  (lib) 

T(n  +  l)  v  ' 

Equation  (11)  is  valid  for  v  e(0.5,  1)  <=>  y2  e(3,  4).  For  a<a,hL  =10  works  for  all  valid  values  of  y2.  For 
&  <  a ,  the  values  of  both  a  and  %  depend  on  the  value  of  the  spectral  exponent  y2:  it  has  been  determined  [19] 
that  the  integral  sum  in  Eq.  (11a)  calculates  the  integral  1(a)  accurately  to  three  significant  figures  for  the  a  and 
hjj  given  in  Table  1  for  nine  values  of  y2  e  (3,  4). 

It  is  straightforward  to  show  that,  when  only  the  first  term  of  Eq.  (lib)  is  used,  Eq.  (7)  becomes 


*  in  t  ' 


W0 


Kl Qh 


i\2(v+l) 


W(Qh), 


(12) 


which  is  the  result  from  first-order  perturbation  theory. 


Table  1  —  Sum-Transition  Parameter  Values  for  Eq.  (11a) 


Y: 2 

3.1 

3.2 

3.3 

3.4 

3.5 

3.6 

3.7 

3.8 

3.9 

a 

0.39 

0.25 

0.20 

0.15 

0.12 

0.08 

0.05 

0.03 

0.025 

hjj 

6545 

610 

155 

93 

62 

60 

62 

67 

54 
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Rational-Function  Approach 

To  obtain  an  expression  for  1(a)  that  extends  its  range  of  validity  to  all  v2  e  [0.2,  0.95],  a  second  numerical 
procedure  has  been  developed  [20-21].  This  approach  uses  expansion  (lib)  with  hL  =40,  and  replaces  expan¬ 
sion  (11a)  by  exploiting  a  rational- function  approximation  of  B(y,  a)  for  a  particular  reference  value  of  a.  This 
transition  value  a  depends  on  y2  via  the  empirical  relationship  [18] 

log10  &  =  -1.7567y2  +4.4926,  (13) 

and  replaces  the  discrete  a’s  of  Table  1. 

Defining  x  =  y2,  the  rational- function  approximation  of  order  [/V,  D\,  with  D  >  A,  is 


B(4x,a) 


1  +  nxx-\ —  +  nNxN 
1  +  di  +  ---  +  dDxD 


(14) 


With  the  numerator  and  denominator  polynomials  factored,  the  right-hand  side  of  Eq.  (14)  becomes  (1  -  xlz{) 
•••  (1  -  xlz,N)l(\  -  x/p :)  •••  (1  -  x/pD)  in  terms  of  the  zeros  zt  and  poles  p-.  Any  poles  that  happen  to  lie  on  the 
positive  real  axis  will  degrade  the  quality  of  the  approximation.  These  are  dealt  with  using  a  splitting  technique 
that  replaces  the  pole  p  >  0  with  a  pair  of  rotated  poles  (pe+lh ,  pe~lh)  and  also  introduces  a  new  zero  z  =  pl(2 
cos  8-1).  Each  such  replacement  increases  both  A  and  D  by  unity,  leaving  D  -  A  unchanged.  The  rotation  angle 
8  is  arbitrary  except  for  the  requirement  8  >  60°  to  prevent  z  >  0.  (We  use  8  =  90°. )  After  the  rotation  of  any 
positive  poles,  the  factored  form  of  Eq.  (14)  can  be  converted  to  the  partial-fraction  form, 
=  X  ._x  Tj  /( x-pj ),  involving  the  modified  poles  and  the  residues  r-.  For  general  a,  the  exponential 
factor  can  be  written  B(y,  a)  =  B(sy,  &)  in  terms  of  the  scale  factor  s  =  ( a/a)1/2v  so  that 
B(y ,  a)  =  X  ._x  Vj  /  ((5y)  -  pj) .  With  this,  Eq.  (8)  reduces  to  a  sum  of  tabulated  integrals,  and  the  final  outcome 
is 


■fpj 


(15) 


where  K0  is  the  zeroth-order  modified  Bessel  function  of  the  second  kind.  In  practice,  one  can  use  [A,  D ]  =  [6, 
9]  in  all  cases  [18]. 

To  apply  these  equations  to  the  ocean  bottom  and  surface,  we  need  to  specify  the  appropriate  boundary 
conditions.  This  is  done  in  the  next  two  subsections. 

Ocean  Bottom  Model 

Wurmser  [13]  has  derived  (3  for  the  general  case  of  a  fluid/elastic-solid  interface.  In  this  case,  the 
following  additional  environmental  parameters  are  required  beyond  c0: 

phot  =  ratio  of  solid  mass  density  to  the  fluid  mass  density 
cp  =  complex  sound  speed  of  the  p- wave  in  the  solid  (compressional  speed  and  attenuation) 
cs  =  complex  sound  speed  of  the  s-wave  in  the  solid  (shear  speed  and  attenuation). 

The  complex  acoustic  wavenumbers  corresponding  to  cp  and  cs  are,  respectively: 


kp  =  2k  f  /  cp  and  ks  =  2k  f  /  c, . 


(16) 
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We  now  define  three  pairs  of  general  functions  that  will  be  used  in  constructing  (3: 


where 


a/iOjip^k^k  ,ks)  =  pbot 


k0  sin  Oj 


4 K  cosz  OjSjiSj  - Pj ) 


^ j  j  ’  Pfotf  ’  ^0  ’  )  Pbot 


k0  sin  Oj 


1- 


2£q  cos2  0: 


(17) 


%j(Oj;k0,k  ,ks)  = 


k0  sin0,  r  .  0-i 


P,  =  f  l -kl  cos2  0,  ,  S,-  =  jk2  -kl  cos2  0  ,- , 
with  0,  =  G/nr  and  02  =  0**. 


(18) 


Finally,  we  note  from  Eqs.  (2)  and  (3)  that 

[kh  <lh\  =  k 0  Cos0mc  COS0«:«/  COS^,. . 


(19) 


With  these  definitions,  we  are  now  ready  to  write  down  the  expression  for  (3: 


P(binc’  bscat ’  Pbi  ’  Pbot’  k0  ’  ks  ’  kp ) 
.  ^Pbot 


p'  (ax+l)(a2+l) 


4(1  -  P*  j([*A  '  <Ih  ]bi4  -  *o2  sin  0/«c  sin0*J 


2  kt 


4k2p 


k ? 


[*/,  •?*]  -4[*A  ?a](«1  -Il)(«2  -|2)+4^o«1«2  +  -^-(«1  _  (’l  )(a2  ~  bl) - ~  blb2 

Pbot  Pbot 


(20) 


Thus,  for  this  model,  depends  on  eight  environmental  parameters:  c0  and  seven  geophysical  param- 
eters:  pbot,  Re{cp},  lm{cp},  Re{cs},  Im{c,},  w2,  and  y2. 

It  should  be  noted  that  this  3D  formulation  of  a- b°f  is  for  a  fluid-solid  interface  (and  reduces  to  the  fluid- 
fluid  case  in  the  appropriate  limits).  However,  for  most  cases  this  SSA  formula  will  handle  not  only  the  water- 
sediment  interface,  but  sediment-sediment  and  sediment-(acoustic)-basement  interfaces  as  well.  For  these  latter 
scenarios,  the  overlying  sediment  basically  acts  as  a  fluid,  so  that  the  above  expressions  would  still  hold  for 
such  cases  as  a  sand-limestone  interface.  (For  expressions  capable  of  handling  solid-solid  interfaces,  see  Ref. 
13.) 

Ocean  Surface  Model 

Calculating  P  for  the  sea  surface  is  much  simpler  (Dirichlet  boundary  conditions).  The  result  is  [12] 


P  =  -4^o  sin0,„c  sin0«.a; . 


(21) 
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(Note  that,  in  contrast  to  the  ocean  bottom  case,  (3  for  the  ocean  surface  does  not  depend  on  the  bistatic  angle 
§bi-) 


The  sea  surface  contains  many  scales  of  roughness,  from  the  long  gravity  waves  (“swell”)  to  the  short 
capillary  waves.  While  a  variety  of  directional  surface-roughness  spectral  models  are  available  [22],  for  this 
paper  we  assume  an  isotropic  spectral  model  (Eq.  (6)),  with  a  spectral  strength  given  by 

w2=AsU,  (22) 

where  U  is  the  wind  speed  (in  m/s)  at  an  elevation  of  10  m.  With  this  W(Qh)  model,  for  a  given  wind  speed  there 
are  two  free  parameters  with  which  to  fit  data  curves,  y2  and  As,  y2  driving  the  frequency  dependence,  and  y2 
and  As  the  level.  Typical  open-ocean  values  of  these  parameters  are:  y2  e(3.4,  4)  and  As  e(5  x  10-5,  20  x 
10-5)  m3-s.  (More  general  forms  of  Eq.  (6)  are  under  consideration  and  will  be  incorporated  into  future 
models.) 

Thus,  for  this  model,  depends  on  four  environmental  parameters:  c0,  U,  As ,  and  y2.  Note,  with  this 
model,  depends  on  frequency  only  through  its  dependence  on  y2,  and  depends  on  the  bistatic  angle  only 
through  its  dependence  on  \Qh\. 

4.  BOTTOM  INTERFACE  SCATTERING 

Our  3D  bottom-interface  scattering  model  was  described  in  the  previous  section.  Here,  after  a 
discussion  of  its  environmental  inputs,  we  briefly  explore  some  of  its  implications. 

Geophysical  Parameter  Values 

Table  2  lists  some  representative  model  parameter  values  for  a  range  of  bottom  types:  mud  (silty  clay), 
(coarse)  sand,  mudstone,  limestone,  and  (rough)  basalt.  (These  values  derive  from  a  number  of  sources.  The 
parameters  for  mud,  sand,  and  basalt  derive  from  Essen  [23]  and  Jackson  [24];  those  for  limestone  derive  from 
Hamilton  [25,26]  and  Jackson  [24];  while  those  for  mudstone  derive  from  Fulford  [27]  and  Kunz  and  Gauss 
[28].) 


Table  2  —  Representative  Bottom-Interface  Seafloor  Parameters 


P  bot 

cP/c 0 

72 

w2  (m4) 

Mud 

1.4 

0.99-0.002/ 

0.1-0.004/ 

3.3 

0.000518 

Sand 

2.0 

1.2-0.005/ 

0.3-0.07/ 

3.3 

0.006957 

Mudstone 

2.2 

1.5-0.001/ 

0.5-0.0005/ 

3.5 

0.0024 

Limestone 

2.4 

2.7-0.004/ 

1.5-0.003/ 

3.3 

0.004 

Basalt 

2.7 

3.4-0.006/ 

1.8-0.006/ 

3.3 

0.01862 

Table  2  illustrates  some  typical  parameter  value  trends,  such  that  as  the  hardness  of  the  bottom  material 
increases,  the  density  and  sound  speeds  increase.  Another  observation  is  that  the  shear  speed  is  often  much  less 
than  the  compressional  speed.  Indeed,  the  ratio  of  shear  speed  to  compressional  speed  is  constrained  as  follows 
[25]: 


Re{c,}  <  1 1-2\|/ 
Refcp}  _  \  2  —  2\|/  ’ 


(23) 
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where  \|/  is  Poisson’s  ratio.  Typically  [25],  \| /  ranges  between  0.2  (hard  materials)  and  0.5  (soft  materials— \|/  = 
0.5  for  water),  so  that  the  corresponding  range  for  this  ratio  is  between  0.6  and  0.  As  regards  bottom  roughness, 
unfortunately  (in  contrast  to  the  sea  surface)  few  measurements  exist.  In  view  of  this  incomplete  knowledge,  it 
is  typically24  assumed  that  y2  ~  3.3,  and  w2  generally  increases  with  material  hardness.  In  practice,  when  in  situ 
acoustic  data  are  available,  these  latter  two  parameters  are  often  estimated  empirically  via  geoacoustic  inversion 
(e.g.,  see  Ref.  29.) 

For  the  model  illustrations  in  this  paper,  the  imaginary  parts  of  the  sound  speeds  (which  induce  attenua¬ 
tion)  will  be  assumed  to  be  independent  of  frequency.  Since  our  computations  of  depend  on  the  acoustic 
wavenumbers  (Eq.  (16))  (not  explicitly  on  the  sound  speeds),  this  implies  that  the  compressional-  and  shear- 
speed  attenuations  will  increase  linearly  with  frequency  [26] .  (It  should  be  noted  that  although  this  is  a  standard 
assumption,  there  is  evidence  in  the  literature  to  the  contrary,  especially  for  water-saturated  sediments  [30-31]. 
This  could  have  a  significant  impact  on  sediment  acoustics,  especially  at  mid-frequencies  [32].)  Conversion 
from  attenuation’s  description  as  the  imaginary  part  of  the  wave  velocities  Im {c;},  where  j  =  p,s ,  to  the  more 
commonly  defined  bottom  sound  absorption  coefficient  a :•  =  k  ■  /(kHz)  is  provided  by  [8] 


4000071  Im  {c^ 
ln(10)  [r  e{Cj}]2 


(24) 


If  a j  is  in  its  standard  units  of  dB/m  and  the  sound  speeds  in  m/s,  then  the  acoustic  frequency  must  be  in  kHz 
(hence,  the  subscript).  (For  example,  for  the  sand  of  Table  2,  as  cp  =  (1.2  -  0.005/)c0,  we  have  Kp  =  189.5/ 
c0,  or  0.125  dB/m/kHz  for  c0  =  1515  m/s.) 

Model  Illustrations 

We  now  illustrate  the  model  by  using  parameters  covering  a  range  of  environments,  frequencies,  and 
geometries.  (For  all  modeling  runs,  we  assumed  c0  =  1515  m/s  and,  except  for  Figs.  4(e)  and  4(f),  §bi  =  180°.) 
While  a  comprehensive  parameter  study  is  beyond  the  scope  of  this  paper,  examples  are  provided  to  illustrate 
both  how  bottom-interface  scattering  strength  depends  on  these  quantities  and  how  complex  the  scattering 
response  can  be.  This,  in  turn,  argues  for  the  use  of  physics-based  models  (in  lieu  of  an  empirical  model,  such 
as  that  of  Ellis  and  Crowe  [33])  as  bottom  scattering  strength  submodels  to  active  performance  models,  such  as 
that  of  Fromm  et  al.  [34],  for  predicting  bistatic  reverberation. 

Model-Model  Comparisons 

Figure  2  illustrates  some  key  points  regarding  the  use  of  both  SSA  and  standard  first-order  perturbation 
theory  (PT)  models.  In  this  figure,  monostatic  predictions  of  backscattering  strength  at  3  kHz  are  plotted  as  a 
function  of  grazing  angle.  For  the  mud  case,  we  see  that  all  three  predictions  (SSA,  PT  for  a  fluid-solid 
interface,  and  PT  for  a  fluid- fluid  interface)  lie  atop  one  another  (except  near  specular).  However,  as  we 
progress  to  harder  materials,  elastic  effects  become  increasingly  important  below  the  interface  critical  angles 
(QJcr  =cos-1(c0 /Re{cy})  when  c0  /  Re{cy}<  1,  for  j  =  p,  s).  (Recall,  Re{c5}  <  0.6  •  Re  {cp}  so  that 
0£r  <  9£.  (for  Qscat  <  90°)  when  both  exist.)  In  these  examples,  for  the  sand  and  mudstone  cases,  fluid-fluid  PT 
generally  overpredicts  the  backscattering  strengths  below  the  compressional-speed  critical  angles  (-33.5°  and 
49°,  respectively),  while  for  both  the  basalt  and  limestone  cases,  fluid-fluid  PT  generally  underpredicts  the 
backscattering  strength  at  low  angles  (below  -40°).  In  contrast,  the  SSA  and  fluid-solid  PT  predictions  are 
relatively  close  below  the  compressional-speed  critical  angles  in  all  cases  except  for  basalt,  where  there  is  only 
agreement  below  -40°.  These  examples  argue  that  an  elastic  model  is  needed  to  accurately  predict  the  acoustics 
at  low  scattering  angles. 

The  basalt,  limestone,  and  mudstone  cases  also  illustrate  the  role  that  elastic  effects  can  play  in  complicat¬ 
ing  the  scattering  response  in  the  vicinity  of  and  above  the  critical  angles.  For  example,  in  the  limestone  case  of 
Fig.  2(b),  the  local  peak  at  49°  and  null  at  68°  correspond  to  the  shear  and  compressional  critical  angles, 
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0  10  20  30  40  50  60  70  80  90 

GRAZING  ANGLE  (deg) 


Fig.  2  —  Monostatic  model  comparisons  of  the  grazing-angle 
dependence  of  backscattering  strength  at  3  kHz  using  the  val¬ 
ues  in  Table  2.  Included  are  the  SSA  and  PT  predictions  for  (a) 
mud,  sand,  and  basalt,  and  (b)  mudstone  and  limestone.  For 
reference,  a  curve  corresponding  to  the  Mackenzie  rule  is  in¬ 
cluded  in  each  plot.  The  p’s  and  s’s  identify  the  compres- 
sional-  and  shear- speed  critical  angles,  respectively. 


respectively.  Furthermore,  when  the  shear  speed  is  significant,  substantial  subcritical-angle  nulls  can  develop. 
This  can  occur  even  when  Re{cJ  is  significantly  less  than  c0  as  demonstrated  with  the  mudstone  case,  where 
the  shear  speed  is  half  the  water  sound  speed.  This  example  reinforces  the  importance  of  using  an  elastic  model. 

Another  feature  of  the  SSA  is  that  its  predictions  behave  well  at  all  scattering  angles,  whereas  PT  begins  to 
have  problems  near  the  specular  direction  (typically  overpredicting  the  scattering  strength)  (Fig.  2).  Also,  note 
that  the  peak  scattering  response  is  not  always  at  specular.  Although  not  illustrated,  rough-surface  scattering 
strength  can  be  quite  sensitive  to  the  value  of  the  spectral  exponent  y2  [18],  so  that  even  for  moderate  grazing 
angles,  the  SSA  is  required  to  capture  the  physics. 

Figure  2  highlights  the  advantage  of  using  the  best  physics-based  models  available,  as  distinguished  from 
the  Navy- standard,  phenomenological  approaches  (i.e.,  a  10  •  log10  (p  sina0^c  sina  Qscat)  dependence  for  a  = 
0.5  or  1,  where  p  is  determined  empirically).  To  illustrate  this,  the  Mackenzie  curve  (-27  +  10  •  log10  sin2  0 
dB)  has  been  included  in  Fig.  2  to  allow  comparison  to  the  SSA  model  predictions.  It  is  seen  that,  even 
translating  the  curves  vertically,  any  agreements  between  the  SSA  and  one-parameter  empirical  predictions  are 
usually  limited  to  a  small  set  of  angles.  Hence,  if  acoustic  data  were  just  available  at  intermediate  angles  (say, 
15-30°)  and  fit  with  an  empirical  curve,  then  one  would  typically  overestimate  the  backscattering  contributions 
at  low,  “long-range”  grazing  angles  (e.g.,  5°).  Another  general  danger  of  the  phenomenological  approach  in 
extrapolating  data  values  to  lower  (or  higher)  angles  is  that  the  scattering  mechanism  can  change  (e.g.,  sedi¬ 
ment-volume  scattering  above  the  compressional  critical  angle,  and  interface  scattering  below). 
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Model  Predictions  for  Five  Bottom  Types 

Figure  2(a)  shows  the  wide  range  of  backscattering  values  (30  dB  or  more  at  a  given  angle)  that  can  be 
encountered  in  the  ocean.  Moreover,  this  variance  depends  strongly  on  the  grazing  angle  (and  frequency).  In 
this  example  at  3  kHz,  Fig.  2(a)  shows  that  at  5°  the  acoustic  response  from  basalt  is  about  35  dB  stronger  than 
that  from  mud  or  sand;  at  20°  the  acoustic  response  from  basalt  is  about  15  dB  stronger  than  that  from  sand, 
which  in  turn  is  20  dB  stronger  than  that  from  mud;  while  at  80°  the  acoustic  responses  are  roughly  equal. 
Although  these  predictions  are  for  3  kHz,  the  general  implications  hold  for  a  wide  range  of  frequencies— e.g., 
model  predictions  at  300  Hz  qualitatively  look  very  similar  to  those  of  Fig.  2,  so  that,  except  toward  specular, 
the  relative  backscatter  strengths  for  the  five  materials  described  by  Table  2  are  generally  maintained  over  10  to 
10000  Hz.  This  is  illustrated  in  Fig.  3,  which  shows  the  frequency  dependence  of  their  backscattering  strengths 
as  predicted  by  the  SSA  at  three  grazing  angles.  This  figure  suggests  that  interface  scattering  strength  generally 
has  a  monotonic  frequency  dependence,  with  levels  typically  increasing  moderately  with  increasing  frequency 
(the  rate  depending  on  the  value  of  y2).  An  exception  to  this  rule  of  the  thumb  occurs  at  high  angles,  where  the 
scattering  strength  can  flatten  and  even  begin  decreasing  as  the  frequency  increases. 


FREQUENCY  (Hz) 


Fig.  3  —  Monostatic  predictions  of  the  frequency 
dependence  of  backscattering  strength  at  three  graz¬ 
ing  angles  for  mud,  sand,  mudstone,  limestone,  and 
basalt  using  the  geophysical  values  in  Table  2. 
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Ecstatically,  similar  trends  to  those  discussed  above  emerge.  This  is  illustrated  in  Fig.  4,  where  the 
dependence  of  scattering  strength  on  the  scattered  grazing  angle  dscat  is  highlighted.  Figure  4(a)  presents  a 
bistatic  version  of  Fig.  2(a)  for  an  incident  grazing  angle  Qinc  of  10°,  suggesting  that  the  general  trends  in 
scattering  strength  both  between  models  and  bottom  types  are  preserved.  (Also  note  that,  bistatically  the  critical 
angles  occur  at  QJcr  and  180°  -QJcr,  where  j  =  p,s.) 


10  50  90  130  170 

SCATTERING  ANGLE  (deg) 
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10  50  90  130  170 
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-40  -31  -22  -14  -5 
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Fig.  4  —  Bistatic  scattering  strength  predictions  for  (a)  mud,  sand,  and  basalt,  and  (b)-(f) 
limestone,  as  a  function  of  dscat  using  the  geophysical  values  in  Table  2.  (The  legend  of  Fig. 
2  applies  to  (a).)  The  black  dashed  lines  in  (a)-(d)  and  the  S’s  in  (e)-(f)  indicate  the  specular 
direction  in  each  case,  while  the  white  lines  in  (b)-(d)  and  the  B’s  in  (e)-(f)  indicate  the 
backscatter  directions.  The  p’s  and  s’s  in  the  margins  of  (b)-(f)  indicate  the  compressional- 
and  shear- speed  critical  angles,  respectively. 
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To  explore  this  assertion  more  systematically,  Figs.  4(b)  through  4(f)  display  SSA  predictions  of  scattering 
strength  as  a  function  of  frequency  and  the  three  scattering  angles  for  a  single  bottom  type  (the  limestone  of 
Table  2).  Figure  4(b)  shows  scattering  strength  at  3  kHz  as  a  function  of  both  the  incident  and  scattered  grazing 
angles.  The  critical-angle  effects  are  readily  apparent,  in  this  case  when  Qinc  or  Qscat  =  -50°  or  -130°  (shear- 
speed  critical  angles),  and  when  Qinc  or  Qscat  =  -68°  or  -112°  (compressional- speed  critical  angles).  As 
expected,  we  see  that  the  scattering  response  is  strong  near  specular.  (We  also  see  the  requisite  symmetry  about 
Vine  =  ®scat  for  ®scat  ^  90°  and  about  einc  =  1  80 °  “  ®scat  for  90°  <  0,raf  <  180°.) 

Figures  4(c)  and  4(d)  examine  the  scattering  response  at  two  incident  grazing  angles  as  a  function  of 
frequency  and  the  scattered  grazing  angle.  We  see  that  for  a  fixed  value  of  QscaP  the  scattering  strength  generally 
increases  with  increasing  frequency  at  low  scattering  angles,  but  not  always  at  high  scattering  angles.  Figures 
4(e)  and  4(f)  examine  the  scattering  response  at  3  kHz  for  the  same  two  incident  grazing  angles  as  a  function  of 
the  bistatic  and  scattered  grazing  angles.  Not  unexpectedly  given  our  isotropic  roughness  assumption,  away 
from  the  vicinities  of  the  specular  direction  and  the  critical  angles,  the  dependence  on  bistatic  angle  is  relatively 
mild.  (We  also  see  the  requisite  symmetry  about  §bi  =  0°,  and  that  SS(Qinc,  Qscat ,  §bi)  =  SS(Qinc,  180°  -  Qscat,  <\>bi 
±  180°).) 

Figure  4  demonstrates  that  bottom  scattering  is  a  complex  process  requiring  physics-based  models  to 
properly  estimate  its  acoustic  response. 

Environmental  Parameter  Dependencies 

Figures  5  and  6  present  a  study  of  the  sensitivity  of  scattering  strength  to  geophysical  parameter  values  as 
functions  of  grazing  angle  and  frequency,  respectively.  The  goal  is  not  to  provide  a  complete  depiction  of  how 
Gfa  depends  on  its  inputs,  but  rather  to  gain  a  feel  for  the  relative  angle-  and  frequency-sensitivity  of  o^f  to 
geophysical  parameters.  The  scheme  adopted  for  this  purpose  was  to  use  the  sand  parameters  of  Table  2  as  the 
default,  and  then  independently  vary  the  seven  geophysical  parameters  over  a  wide  range  of  realistic  values.  For 
simplicity,  the  study  was  restricted  to  monostatic  SSA  predictions  at  1500  Hz. 

Figure  5  reveals  both  how  interesting  the  grazing-angle  behavior  of  off/  can  be,  and  its  sensitivity  to  its 
input  geophysical  parameter  values.  It  is  apparent  that  across  most  of  the  figures,  scattering  strength  generally 
increases  as  the  grazing  angle  increases.  Exceptions  to  this  trend  are  regularly  seen  in  the  vicinities  of  the 
critical  angles.  For  example,  Fig.  5(a)  shows  a  two-peaked  ridge  of  maximum  values,  running  from  the  vicinity 
of  the  critical  angle  at  high  w2  values  to  specular  at  low  w2  values.  Figure  5(b)  shows  that,  for  a  given  grazing 
angle,  scattering  strength  has  a  fairly  flat  response  as  a  function  of  compressional  speed,  above  the  ridge 
tracking  the  compressional-speed  critical  angle  (the  value  of  which  increases  as  the  compressional  speed 
increases).  Below  this  ridge,  this  situation  is  different:  for  a  given  compressional  speed,  the  scattering  strength 
has  a  mild  monostatically-increasing  response  as  a  function  of  grazing  angle. 

Figures  5(c)  and  5(d)  reveal  grazing  angle  valleys  for  particular  combinations  of  shear  speed  and  grazing 
angle  values.  To  provide  a  more  complete  picture  of  the  role  of  shear  in  influencing  the  acoustics,  the  default 
compressional  speed  of  1818  m/s  has  been  increased  to  4500  m/s  in  Fig.  5(d).  In  this  case,  a  shear- speed 
critical  angle  exists  when  Re{cJ  >  c0  (as  mapped  out  by  the  white  line).  As  a  result,  a  second  valley  tracking 
the  shear-speed  critical  angle  has  emerged.  As  illustrated  earlier  with  the  mudstone  case  (Fig.  2(b)),  these 
shear-induced  nulls  are  significant  manifestations  of  an  elastic  theory  that  needs  to  be  accounted  for  when 
modeling  reverberation  in  rocky  environments  (see,  e.g.,  Ref.  28).  Whether  appreciable  nulls  manifest  them¬ 
selves  for  a  given  bottom  at  a  given  frequency  will  depend  on  the  bottom’s  geophysical  values.  For  example,  in 
our  artificial  bottom  of  Fig.  5(d),  if  the  shear  speed  were  1800  m/s,  there  would  be  a  significantly  weaker 
backscattering  response  at  low  grazing  angles  than  if  the  shear  speed  were  2300  m/s. 

For  a  given  grazing  angle,  Figs.  5(e),  5(g),  and  5(h)  suggest  a  relatively  mild  dependence  of  scattering 
strength  on  the  roughness  spectral  exponent  and  the  attenuation  values,  while  Fig.  5(f)  shows  the  expected 
increase  in  scattering  strength  as  the  density  ratio  increases.  Note  in  these  four  cases,  local  maxima  in  the 
vicinity  of  the  compressional  critical  angle,  and  that,  except  for  the  density  ratio,  these  maxima  are  typically 
strongest  for  low  parameter  values,  with  these  maxima  getting  “smoothed  out”  for  sufficiently  high  parameter 
values.  For  more  elastic  bottoms,  a  stronger  dependence  on  y2  would  be  observed  [18]. 
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Fig.  5  —  Monostatic  predictions  at  1500  Hz  of  the  grazing  angle  dependence  of  bottom  backscattering 
strength  on  geophysical  parameter  values.  Using  the  sand  parameters  of  Table  2  as  the  default,  the  displayed 
geophysical  parameters  were  independently  varied  over  the  ranges  shown.  In  (d),  the  default  value  of 
Re{cp}  was  reset  to  4500  m/s,  before  varying  the  shear  speed.  The  p’s  and  the  black  curve/arrow  in  (b) 
indicate  the  compressional-speed  critical  angles  in  each  case,  while  the  s  and  the  white  curve  in  (d)  indicate 
the  shear-speed  critical  angles.  The  gray,  dashed  horizontal  lines  demark  the  default  parameter  values. 
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Fig.  5  (continued)  —  Monostatic  predictions  at  1500  Hz  of  the  grazing  angle  dependence  of  bottom  backscat- 
tering  strength  on  geophysical  parameter  values.  Using  the  sand  parameters  of  Table  2  as  the  default,  the 
displayed  geophysical  parameters  were  independently  varied  over  the  ranges  shown.  The  p’s  indicate  the  com- 
pressional-speed  critical  angle  in  each  case.  The  gray,  dashed  horizontal  lines  demark  the  default  parameter 
values. 
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Fig.  6  —  Monostatic  predictions  of  the  frequency  dependence  of  bottom  backscattering  strength  on 
three  geophysical  parameter  values  at  two  grazing  angles.  Using  the  sand  parameters  of  Table  2  as  the 
default,  the  displayed  geophysical  parameters  were  independently  varied  over  the  ranges  shown.  The 
gray,  dashed  horizontal  lines  demark  the  default  parameter  values. 
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Figure  6  shows  some  representative  trends  in  the  frequency  behavior  of  bottom  interface  scattering  strengths 
(cf.,  Fig.  3).  Most  geophysical  parameters  influence  G^f  in  the  manner  of  the  density  ratio,  namely,  a 
monotonically  increasing  dependence  on  frequency  for  a  fixed  value  of  the  parameter  at  low-to-moderate 
grazing  angles,  with  a  more  complex  behavior  at  higher  angles,  for  example,  the  ridges  of  maximum  scattering 
strengths  at  60°  seen  in  Figs.  6(e)  and  6(f).  Figures  6(c)  and  6(f)  display  the  expectedly  weaker  frequency 
dependence  as  the  value  of  the  spectral  exponent  gets  closer  to  4. 

Although  not  intended  to  be  a  comprehensive  look  at  the  dependence  of  of^f  on  its  input  parameters,  this 
parameter  study  does  suggest  that,  outside  of  the  vicinity  of  the  critical  angles,  the  spectral  parameters,  sound 
speeds  and  density  will  generally  have  the  most  influence  on  the  acoustics  at  low  grazing  angles  (at  least  for  y2 
e  [3,  4]).  Hence,  obtaining  good  estimates  for  their  values  will  be  especially  important. 

Model-Data  Comparison 

Figure  7  shows  an  example  data-model  comparison  of  for  a  rocky  bottom.  The  monostatic  data  were 
collected  over  an  exposed  limestone  bottom  off  the  Carolina  coast  during  the  LWAD  Focused  Technology 
Experiment  (FTE)  96-2  [35-36].  In  this  case,  below  the  compressional  critical  angle,  the  model  matched  the 
observed  grazing  angle  dependence,  obtaining  an  rms  data/model  mismatch  of  1.3  dB  over  backscatter  angles 
of  15  to  70°  for  four  frequencies  across  the  2.0  to  3.5  kHz  band  [36].  The  unknown  bottom  parameters  were 
estimated  using  simulated  annealing:  pbot  =  2.6,  cp  =  (3.2  -  0.12/)c0,  cs  =  (1.5  -  0.48z)c0,  w2  =  0.0004  m4, 
and  y2  =  2.6.  (In  this  case,  c0  =  1487  m/s.) 

This  example  illustrates  how  the  model  can  be  used  to  invert  for  hard-to-measure  roughness  spectral 
parameters,  in  particular  for  y2:  since  it  is  the  only  environmental  parameter  that  has  any  impact  on  the 
frequency  dependence,  broadband  measurements  can  be  used  to  estimate  it.  (Note,  this  assumes  that  the 
complex  sound  speeds  can  be  taken  to  be  effectively  independent  of  frequency  over  the  measurement  band.) 
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Fig.  7  —  Data-model  comparison  for  an  exposed 
limestone  bottom  off  the  South  Carolina  coast.  The 
shallow-water  acoustic  data  (symbols)  were  collected 
monostatically  using  an  omnidirectional  source  and 
a  vertical-line  array  receiver  during  LWAD  FTE  96- 
2  at  Site  Q  (from  Ref.  36). 
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In  comparing  these  curves  with  the  limestone  curve  in  Fig.  2(b),  note  how  in  these  LWAD  96-2  model 
predictions,  the  increased  attenuation  values  have  smoothed  out  the  local  peaks  and  nulls  associated  with  the 
critical  angles.  Note  also  the  model  underestimates  the  data  above  the  compressional  critical  angle,  where 
volume-scattering  contributions  are  expected. 

Comments 

While  Figs.  2  through  7  have  revealed  complex  relationships  between  and  its  input  parameters,  these 
can  be  used  to  advantage  in  inverting  for  geophysical  parameter  values.  By  maximizing  the  frequency  and 
grazing-angle  (both  incident  and  scattered  if  possible)  coverage,  systematic  identification  of  these  values  can  be 
pursued  by  exploiting  distinctive  acoustic  signatures  (such  as  critical-angle  and  shear- speed  effects).  Typically 
for  a  given  ocean  bottom,  some  values  (such  as  the  density  ratio)  are  known  with  more  confidence  than  others 
(such  as  the  spectral  strength).  Hence,  the  use  of  a  physics-based  model  such  as  the  one  presented  here  will 
allow  an  improved  estimation  of  these  other  values,  some  of  which  can  be  difficult  to  measure  even  with  the 
most  sophisticated  equipment.  Furthermore,  if  these  geophysical  parameters  are  estimated  from  a  small  band  of 
angles  at  a  few  frequencies,  use  of  such  a  model  will  allow  extrapolation  in  frequency  and  scattering  angles  with 
a  degree  of  confidence  that  is  not  possible  with  empirical  (e.g.,  Lambert’s  law)  approaches.  A  physics-based 
acoustic  model  can  also  provide  confidence  levels  to  its  predictions  based  on  the  confidence  levels  associated 
with  the  knowledge  of  its  environmental  inputs. 

The  model  presented  describes  scattering  from  the  water-sediment  interface  (as  well  as  from  sediment- 
sediment  and  sediment-basement  interfaces).  Not  described  in  this  paper  is  a  general  model  for  sediment- 
volume  scattering.  For  such  models,  the  reader  is  referred  to  the  work  of  Jackson  [8,16,37-40].  Coupling  such 
a  model  with  the  interface-scattering  model  presented  here  would  provide  a  general  description  of  bottom 
scattering,  which  in  turn,  could  be  used  to  identify  under  what  conditions  the  dominant  scattering  mechanism 
is  the  interface  or  the  sediment  volume. 

5.  VOLUME  SCATTERING 

Because  of  their  variety  and  dynamic  nature,  estimating  the  scattering  contributions  of  fish  is  particularly 
challenging.  When  fish  are  well  separated  from  the  ocean  surface  or  bottom,  recognizable  broadband  acoustic 
signatures  identifying  their  presence  and  strength  have  been  observed  [41-44].  However,  when  fish  are  in  the 
vicinity  of  an  ocean  boundary,  as  is  common  in  the  littoral,  these  characteristic  fish  signatures  can  undergo 
significant  modification  due  to  boundary-interference  effects  [5,45-46]. 

We  now  describe  physics-based  volume-scattering  models  that  handle  these  boundary-interference  effects, 
and  provide  a  new  understanding  of  the  relative  significance  of  surface,  fish,  and  bottom  scattering  as  functions 
of  frequency,  scattering  angle,  and  biological  and  environmental  factors.  These  plane-wave  models  are  based  on 
an  established  model  of  Love  [2].  We  begin  with  a  review  of  his  target  strength  formula,  and  follow  with  the 
extensions  applicable  to  dispersed  bladdered  fish  near  the  ocean  surface  and  ocean  bottom. 

Single  Fish  Swimbladder  Model 

At  low  to  mid  frequencies,  the  primary  scattering  mechanism  of  a  fish  is  its  swimbladder  [2].  Here,  we 
first  describe  some  of  its  physical  properties,  and  then  describe  its  free-field  acoustic  scattering  properties. 

Swimbladder  Properties 

Swimbladders  resemble  prolate  spheroids  or  cylinders  with  major-to-minor  axis  (length-to-diameter)  ratios 
up  to  10.  However,  near  resonance,  it  has  been  shown  that  the  acoustic  response,  whether  assuming  a  spherical, 
spheroidal,  or  cylindrical  air-bubble  model,  generally  varies  by  at  most  30%  [47-48].  Accordingly,  as  a  first 
approximation  we  consider  a  spherical  swimbladder,  and  then  introduce  corrections  to  account  for  real  fish. 
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The  fundamental  (monopole)  resonance  frequency  fR  of  a  spherical  swimbladder  is  given  by  [2] 


f2  _  3y aP_ 

JR  .  9  2  ’ 

4nrzptissue 


(25) 


where  ya  is  the  ratio  of  specific  heats  of  air  (~1.4),  P  is  the  ambient  pressure  in  Pa,  r  is  the  swimbladder  radius, 
and  ptissue  is  the  density  of  the  surrounding  tissue  (typically  [2]  taken  as  1050  kg/m3). 

As  a  next  step,  we  need  to  be  more  realistic  on  the  following  points  before  the  above  expression  will 
describe  the  acoustic  response  of  oceanic  fish:  (a)  P ,  of  course,  depends  on  the  (fish)  depth  z  >  0;  (b)  r  does 
too,  but  in  a  species-specific  way;  (c)  actual  swimbladders  are  nonspherical,  and  this  too  is  species-specific.  We 
discuss  these  in  turn. 

(a)  It  is  elementary  that  the  ambient  pressure  has  the  linear  dependence  P  =  A  +  Rz,  where  A  is  the 
atmospheric  pressure  (nominally  1  atm,  i.e.,  101325  Pa),  and  R  is  the  rate  at  which  the  hydrostatic  pressure 
increases  with  depth  (approximately  ^4/10  per  meter).  Thus  we  have 


P 


= 101325 • 


1  +  - 


9.8066 


(26) 


(b)  For  an  inanimate  balloon,  there  would  be  a  simple  dependence  of  r  on  z:  an  isothermal  descent  would 
just  compress  the  balloon  according  to  Eq.  (26)  and  Boyle’s  Law  [z  =>  P  =>  V  =>  r\.  Fish,  however,  control 
compression  of  their  swimbladders  in  a  manner  empirically  described  by  [42] 


^  r0  ^  <  Zcomp 

J  9.8066+  zc  V/3  ^  (27) 

v  9.8066+  z  J  r°  z~z™p’ 

where  r0  is  the  swimbladder  radius  at  the  ocean  surface  and  zcomp  is  the  depth  where  compression  sets  in.  Both 
of  these  parameters  are  species-dependent,  (r  is  now  the  radius  of  the  equivalent  sphere  of  the  same  volume  V.) 
While  it  may  seem  that  we  are  relying  too  heavily  on  the  concept  of  a  spherical  swimbladder,  namely  that  we 
seem  to  be  correcting  the  radius  of  a  supposedly  spherical  swimbladder  based  on  the  behavior  of  actual  fish, 
when  they  in  fact  have  nonspherical  swimbladders.  However,  Eq.  (27)  relates  any  linear  measure  r  of  a 
swimbladder  to  its  surface  value  r0.  The  shape  really  need  not  be  spherical.  But,  when  it  is  not,  Eq.  (25)  needs 
to  be  modified. 

(c)  With  swimbladders  modeled  as  prolate  spheriods,  the  modification  takes  the  form  [47] 


/*=%.  (28) 

where  Weston’s  (depth-  and  species-dependent)  correction  factor  £  is  a  function  of  the  swimbladder’s  minor-to- 
major  axis  ratio  8  (<  1)  given  by  [2] 


C  =  e"iy3  12? /In 


f\ V 

1  -q 


(29) 


with  q  =  y  1-82.  Values  of  8  generally  range  between  0.1  and  0.5,  which  correspond  to  £  values  of  1.25  to 
~1.0.  (£  is  monotonically  decreasing  over  this  range  of  8.)  Note  that  as  8  — »  1,  we  have  £  — »  1,  as  expected. 
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In  most  fish,  bladder  length  is  constant,  held  rigid  in  the  anterior-posterior  direction  and  on  its  dorsal 
surface  by  connective  tissue  and  the  skeleton,  while  bladder  girth  is  free  to  compress  inward  from  the  lateral 
and  ventral  surfaces  that  adjoin  the  pliant  abdominal  cavity.  Hence,  we  assign  fish  a  spheroid-shaped  bladder  of 
constant  length  (usually  LI 3,  where  L  is  the  length  of  the  fish).  For  those  fish  where  the  bladder  compresses  (or 

expands)  during  vertical  movements,  we  determine  its  new  girth  using  asph  =  ^3Vcomp  /  4nbsph ,  where  bsph  and 
asph  are  the  major  and  minor  axis  radii,  respectively,  and  Vcomp  is  the  compressed  volume  calculated  from  the 
r  in  Eq.  (27).  The  new  major  and  minor  axes  are  used  in  the  calculation  of  the  spheroid-shape  correction  (Eq. 
(29)),  while  the  equivalent  radius  of  a  sphere  of  the  same  volume  is  used  in  determining  the  resonance 
frequency  (Eq.  (25)). 

In  summary,  Eq.  (28)  can  be  used  as  the  corrected  form  of  the  fundamental  frequency  (Eq.  (25)),  with  Eq. 
(26)  for  the  pressure,  Eq.  (27)  for  the  radius,  and  Eq.  (29)  for  the  shape-correction  factor.  Note  that  fRe  depends 
primarily  on  the  ambient  pressure  (/|  °c  to  P)  and  the  swimbladder  radius  (/|  inversely  °c  to  r2)  that,  in  turn, 
depend  primarily  on  fish  species,  size,  and  depth. 

Scattering  Model 

Our  scattering  model  will  be  that  of  Love  [2] ,  which  represents  an  individual  swimbladder-bearing  fish  as 
an  air-filled,  viscous  elastic  spherical  shell  with  Weston’s  correction.  This  allows  the  acoustic  cross  section 
6  bladder  (in  m2)  in  the  backscattered  direction  for  a  single  swimbladder  at  the  acoustic  frequency /to  be  written 
in  the  form  [2] 


}  bladder  ' 


fl 

f2H2 


+ 


fBE 

f 


(30) 


where  H  is  a  damping  factor: 


1  =  2n tf2  |  % 

H  f. Reco  nr2fRe ptissue  ’ 


(31) 


where  is  the  viscosity  of  fish  flesh  in  Pa-s,  typically  [2]  taken  to  have  a  value  of  50  Pa-s.  (Love’s  model  [2] 
includes  an  additional  term  in  Eq.  (25)  that  accounts  for  the  effects  of  swimbladder  wall  tension  on  fR  and  a 
thermal  damping  term  in  Eq.  (31).  However,  he  has  shown  [2]  that  these  terms  are  not  generally  of  significance, 
so  they  have  been  omitted  here.) 

The  model  implies  that  a  bladdered  fish’s  target  strength  will  depend  primarily  on  its  size  and  depth. 
Figures  8  and  9(a)  provide  representative  illustrations  of  how  the  free-field  frequency  target- strength  response  of 
a  fish,  in  this  case,  the  salmon,  changes  with  size  and  depth,  respectively.  Figure  8  shows  how  the  salmon 
resonance  pattern  varies  as  its  length  is  reduced  from  0.5  m  to  -0.03  m.  (Here,  we  use  a  further  approximation 
that  provides  an  estimate  of  salmon  swimbladder  size  given  its  length.)  We  see  that  as  fish  size  is  reduced,  its 
target  strength  decreases  and  its  resonance  peak  broadens  and  shifts  upward  in  frequency.  Figure  9(a)  shows  the 
salmon’s  resonance  frequency  increasing  with  depth  due  to  the  ambient  pressure  on  its  swimbladder  increasing 
with  depth  (Eqs.  (25)-(26)).  (Also  note  that  at  a  depth  of  10  m,  there  is  a  shift  in  the  frequency  response  as  the 
salmon’s  swimbladder  begins  compression  (Eq.  (27)),  as  for  this  species  [42],  we  assumed  zcomp  =  10  m.) 
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Fig.  8  —  Predicted  frequency  dependence  of  the  target  strength 
of  a  single,  free-field  salmon  at  a  depth  of  10  m  as  parameterized 
by  swimbladder  size  at  the  surface  (~4%  of  salmon  length)  (£ 
was  taken  to  be  1.12  and  zcomp  to  be  10  m) 
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Fig.  9  —  Predicted  frequency  and  depth  dependence  of  near-surface 
monostatic  backscattering  for  a  single  salmon:  (a)  free-field  target 
strength,  and  at  a  grazing  angle  of  10  deg,  the  (b)  Lloyd’s  mirror 
interference  pattern,  and  (c)  near-surface  scattering  response.  ( zcomp , 
r0,  and  £  were  taken  to  be  10  m,  0.0144  m,  and  1.12,  respectively) 
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Layer  of  Fish  Swimbladder  Model 

The  total  scattering  from  a  layer  of  dispersed  fish  is  simply  the  incoherent  sum  of  scattering  from  the 
individuals.  For  the  modeling  in  this  paper,  we  make  a  further  approximation  by  assuming  a  mean  fish  target 
strength  6bladder  at  each  depth  (calculating  a  weighted  average  of  a  Gaussian  distribution  of  Nr  fish  sizes): 


(32) 


r  i= 1 


where  the  wt  are  normalized  weights.  (In  practice,  taking  Nr  =  5  is  often  sufficient.)  Hence  for  a  given  species, 
the  volume  scattering  strength  of  a  layer  of  dispersed  bladdered  fish  will  depend  primarily  on  their  sizes,  depths, 
and  number  densities  (p ylsh). 

Fish  Near  Boundaries 

Some  fish  species  (such  as  salmon  [42])  are  typically  close  to  the  sea  surface  at  night,  and  away  from  the 
surface  during  the  day;  some  species  are  near  the  surface  day  and  night  (such  as  flying  fish  [49]);  and  some 
species  are  close  to  the  bottom  during  the  day,  and  away  from  the  bottom  at  night  (such  as  some  species  of 
rockfish  [50]).  When  bladdered  fish  are  near  ocean  boundaries,  the  scattering  picture  increases  in  complexity. 
Besides  scattering  from  the  rough  interface  itself,  a  fish  scatters  energy  to  a  receiver  along  multiple  paths  (Fig. 
10).  The  relative  time  delay  of  these  four  paths  generates  a  pattern  of  constructive  and  destructive  interference, 
the  intensity  of  which  depends  strongly  on  the  grazing  angles,  the  distance  of  the  scatterer  from  the  boundary, 
the  acoustic  frequency,  and  for  the  ocean  bottom,  its  physical  properties.  This  can  significantly  alter  the  fish’s 
scattered  intensity  (by  a  factor  between  0  and  16). 
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Fig.  10  —  Primary  near-boundary  scattering 
paths,  as  illustrated  for  generic  (a)  monostatic 
and  (b)  bistiatic  surface  geometries.  (Monostatic 
paths  A-B-C-A  and  A-C-B-A  add  coherently) 
(dashed  fishes  are  image  sources) 
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To  model  scattering  from  dispersed  fish  in  the  presence  of  such  boundaries,  we  exploit  the  fact  that  the 
swimbladder  typically  occupies  ~4  to  5%  of  a  fish’s  volume  [51],  allowing  us  to  model  dispersed  fish  as  point 
scatterers.  Furthermore,  our  models  assume  the  ocean  boundaries  are  flat.  (However,  at  our  frequencies  bias 
effects  due  to  a  rough  surface  can  be  shown  not  to  be  a  significant  problem.)  Coupled  with  our  point-scattering 
assumption,  this  implies  that,  to  first  order,  the  scattering  cross  sections  of  bladdered  fish  near  a  boundary 
depend  on  Qinc  and  QscaP  but  not  §bi.  Finally,  these  models  assume  that  the  fish  layer  lies  between  the  sonar  and 
the  boundary,  and  that  the  sonar  is  in  the  far  field  of  the  fish. 

Fish  Near  the  Ocean  Surface 

Because  the  density  of  air  is  so  small,  the  air-sea  interface  basically  reflects  an  incoming  acoustic  wave  so 
that  its  phase  is  flipped.  The  images  of  sources  and  scatterers  near  the  interface  will  also  be  180°  out  of  phase 
from  the  original.  As  a  result,  ensonification  of  a  point  scatterer  at  depth  z  in  the  vicinity  of  the  sea  surface  (Fig. 
10)  creates  an  acoustic  dipole,  enhancing  its  backscatter  cross  section  by  16  sin4  ( k0z  sin  0)  in  the  monostatic 
case  [52].  Figure  9(b)  provides  a  monostatic  example  of  the  resultant  (Lloyd’s  mirror)  interference  pattern  for 
a  grazing  angle  of  10°. 

To  calculate  the  total  scattered  intensity  due  to  a  layer  of  bladdered  fish  near  the  sea  surface,  we  need  to 
convolve  this  Lloyd’s  mirror  effect  with  the  scattering  due  to  the  fish.  Doing  so  for  bistatic  geometries  yields  for 
fish  in  layer(s)  covering  depths  z1  to  z2  (z2  >  zf): 


Z2  _ 

Gfifh  =  16 {  P fisiM)Ohladder(f, z) sin2  (k0 Z sine,,,,.) sin2  (k0z sinQ,cat )dz .  (33) 

*1 


Thus,  for  this  model,  for  fish  uniformly  distributed  (over  a  given  layer)  depends  on  eight  environmen¬ 
tal  and  biological  parameters:  c0,  z1?  z2 ,  zcomp,  pfish,  r0,  8,  and  Nr. 

A  look  at  Eq.  (33)  suggests  that  when  &ozsin0  is  small,  such  as  when  the  scatterer  is  within  a  few 
wavelengths  of  the  surface,  or  the  grazing  angle  is  very  low,  then  the  Lloyd’s  mirror  pattern  can  modify  the  free- 
field  acoustic  response  by  a  factor  between  0  and  16.  Figure  9(c)  shows  a  monostatic  example  of  this  modifica¬ 
tion  at  0  =  10°  for  the  single  salmon  of  Fig.  9(a).  (Note  the  enhanced  scale  of  Fig.  9(c).)  In  this  case,  very  near 
the  surface  the  salmon  resonance  response  is  suppressed  (minimum  of  the  interference  pattern).  At  deeper 
depths,  the  maxima  of  the  Lloyd’s  mirror  pattern  are  more  closely  spaced  in  frequency  and  depth,  so  that  a  layer 
of  fish  over,  say,  20  to  30  m,  ends  up  mimicking  the  free-field  frequency  response  (at  this  grazing  angle),  albeit 
at  enhanced  levels. 

Fish  Near  the  Ocean  Bottom 


In  contrast  to  the  air-sea  interface,  which  effectively  repels  acoustic  energy,  the  water-sediment  interface 
has  some  rigidity  and  attracts  acoustic  energy.  As  a  result,  whereas  the  case  of  a  point  scatterer  near  the  sea 
surface  leads  to  a  relatively  simple  and  compact  formula,  the  general  formulation  for  the  ocean  bottom  will  be 
more  complicated.  (In  the  presence  of  a  non-pressure-release  boundary,  a  point  scatterer  is  not  a  dipole  but  is 
in  general  described  by  a  series  of  higher-order  moments.)  As  with  the  surface  model,  this  bottom  model  has 
as  its  fish-scattering  kernel  the  cross-section  scattering  formula  of  Love  (Eqs.  (25)-(32)).  However,  for  the 
bottom,  we  assume  a  (flat)  fluid/elastic-solid  interface. 


We  first  define  the  complex  bottom-reflection  coefficients  [13] 


a:  -1 

R,  =-L - 

1  a  :  +  1 


(34) 


where  the  «•  are  given  by  Eq.  (17).  Hence,  ab^h  will  depend  on  phlH,  cp,  cs,  and  c0,  but  not  on  w2  or  y2. 
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With  this  definition,  we  can  now  write  down  our  general  bistatic  scattering  cross  section  formula  for  a  layer 
of  bladdered  fish  dispersed  over  depths  z1  to  z2  in  a  waveguide  of  depth  D: 


G%‘h  =  I  P fisiSVGbiaddeM'- ")|]  +  Ruv:  exp[2ik()  (D  -  z) sin  e,„f. ]|2  |l  +  Rscat  exp[2z'&0  (D-z) sin Qscat l|2 dz .  (35) 

*1 


Thus,  for  this  model,  o^h  for  fish  uniformly  distributed  (over  a  given  layer)  depends  on  fourteen  environ¬ 
mental  and  biological  parameters:  c0,  z1?  z2>  zcomp,  A  P fish’  ro>  Nr,  Re{cp},  Im \cp},  Re{cJ,  Im {cs},  and 

P  bot’ 

To  gain  some  insight  as  to  the  effect  of  the  ocean  bottom  on  ab^h ,  consider  the  limiting  case  of  a  perfectly 
reflecting  bottom  (Neumann  boundary  conditions),  where  Rj  —>  1.  In  this  special  case  of  an  ideally  hard,  flat 
bottom,  Eq.  (35)  reduces  to  the  case  of  a  double  monopole: 


Z2  _ 

Ofi*  =16J  P fish (z)obladder(f, z)  cos2[k0  (D -  z)  sin0;„c ]  cos2[fc0  (D  -  Z) Sin evra, ]dz .  (36) 

Z\ 


This  formula  suggests  that  the  acoustic  scattering  response  of  a  fish  near  the  ocean  bottom  at  low  scattering 
angles  is  generally  less  variable  than  that  of  a  fish  near  the  sea  surface  (~cos4  vs  -sin4  for  the  pure  backscatter 
case).  Furthermore,  in  contrast  to  the  sea-surface  case,  a  fish  on  the  bottom  is  in  a  maximum  of  the  interference 
pattern  (source  and  its  image  are  in  phase).  (For  a  realistic,  nonideal  bottom,  this  maximum  will  generally  be 
slightly  above  the  bottom.)  Hence,  fish  resting  on  the  bottom  not  only  still  can  be  appreciable  scatterers,  but 
potentially  mistaken  for  bottom  scatterers  (e.g.,  surface  features  or  sediment- volume  heterogeneities). 

Figure  11  is  an  ocean-bottom  analog  to  Fig.  9.  In  this  example,  we  are  modeling  the  frequency  and  depth 
dependence  of  the  monostatic  backscattering  response  from  a  single  salmon  in  100  m  of  water  at  separations  of 
0  to  30  m  from  a  sandy  bottom  (Table  2).  Figure  11(a)  displays  the  salmon’s  free-field  target  strength  response, 
while  Fig.  11(b)  displays  the  interference  pattern  at  a  grazing  angle  of  10°.  It  is  seen  that  in  this  case,  a  fish  very 
near  the  bottom  will  generally  be  in  a  maximum  of  the  interference  pattern.  We  also  see  that,  due  to  the 
salmon’s  higher  resonance  frequency  at  these  depths,  the  free-field  frequency  response  is  more  closely  mim¬ 
icked  in  this  near-bottom  case  than  in  the  near- surface  case  (Fig.  9(c))  (as  these  maxima  are  more  closely 
spaced  in  depth  at  the  higher  frequencies).  As  with  the  surface  case,  an  important  boundary-interference  effect 
is  the  significant  enhancement  of  the  free-field  backscattering  strengths.  (Note  the  enhanced  scale  of  Fig.  11(c).) 

Comments 

For  practical  applications  of  this  model,  for  example,  as  kernels  in  reverberation  or  performance  models, 
a  field  of  discrete  fish  is  treated  as  a  spatially  (and  temporally)  uniform  layer,  which  in  turn  is  treated  as  an 
effective  modifier  of  boundary  properties.  Thus,  they  can  be  used  like  a  surface  or  bottom  scattering  strength 
with  no  need  to  introduce  a  separate  layer  into  the  modeled  waveguide. 

Model  Illustrations 

We  now  illustrate  the  models  with  a  few  examples.  Here,  for  simplicity,  we  assume  layers  of  fish  that  are 
uniformly  distributed  in  depth  (and,  locally,  uniformly  distributed  in  range).  (For  all  examples,  we  take  §bi  = 
180°.) 
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Fig.  11  —  Predicted  frequency  and  depth  dependence  of  near-bot¬ 
tom  monostatic  backscattering  for  a  single  salmon  above  a  sandy 
bottom  in  100  m  of  water:  (a)  Free-field  target  strength,  and  at  a 
grazing  angle  of  10  deg,  the  (b)  interference  pattern  and  (c)  near¬ 
bottom  scattering  response.  (. zcomp ,  r0,  and  £  were  taken  to  be  10  m, 
0.0144  m,  and  1.12,  respectively.) 
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Fish  Near  the  Ocean  Surface 

The  model  for  ,  when  combined  with  a  surface- scattering  model,  provides  a  prediction  of  scattering 
strength  appropriate  for  the  air-sea  interaction  zone  (when  the  acoustic  source  and  receiver  are  in  the  far  field 
of  the  scatterers).  The  total  scattering  strength  is  the  sum  of  contributions  from  the  air-sea  interface  fish 

(ojsf).  and  tabbies  (oj: 


=  101og10(a"/+c5J+a„).  (37) 

(A  model  for  abub  will  be  described  in  the  section  on  surface  scattering.  Here,  we  take  ahuh  =  0.) 

In  analyses  of  acoustic  backscatter  collected  during  various  low-frequency  (70  to  1500  Hz)  CST  surface- 
scatter  measurements,  the  presence  of  fish  has  sometimes  been  inferred,  particularly  during  periods  of  low  sea 
states  [45,53].  In  these  cases,  for  fish  deeper  than  ~10  m,  recognizable  and  robust  broadband  volume  scattering 
signatures  were  typically  seen  at  low  grazing  angles,  such  as  a  frequency  resonance  and  a  flatness  in  grazing 
angle.  When  fish  are  shallower,  these  acoustic  signatures  can  undergo  significant  modification.  For  example,  in 
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some  of  the  lowest  sea  state  cases  from  the  CST  Gulf  of  Alaska  experiments  (CST-4  and  CST-7  Phase  2), 
measured  low-grazing-angle  scattering  strengths  were  suggestive  of  air- water  interface  scattering,  yet  concur¬ 
rently  measured  continuous  wave  (cw)  spectra  showed  Gaussian  spreads  centered  at  zero  frequency  shift, 
consistent  with  fish  or  bubbles  [45].  However,  no  significant  inclusion  of  bubbles  was  observed  (or,  in  general, 
expected  for  wind  speeds  less  than  5  m/s)  [54] .  The  model  for  just  described  helps  resolve  some  of  these 
mysteries  and  provides  a  new  understanding  of  the  relative  significance  of  interface,  fish,  and  bubble  scattering 
as  a  function  of  frequency,  scattering  angle,  and  biological  and  environmental  factors. 

Figure  12  presents  monostatic  backscatter  model  predictions  for  salmon  uniformly  distributed  over  two  sets 
of  layer  depths,  corresponding  to  their  typical  nighttime  (0.5  to  2  m)  and  daytime  (15  to  40  m)  depths  during 
CST-7  Phase  2  [42].  (Typical  Gulf  of  Alaska  sizes  and  densities  were  assumed  [42].)  Figure  12(b)  displays  the 
frequency  resonance  and  flat  grazing  angle  response  that  occur  when  the  salmon  are  a  number  of  wavelengths 
away  from  the  sea  surface.  Adding  on  the  scattering  contribution  of  the  air-sea  interface,  Fig.  12(d)  reproduces 
the  aforementioned  acoustic  behavior  observed  in  the  CST  Gulf  of  Alaska  surface  scatter  measurements.  In 
contrast,  when  the  fish  are  very  near  the  air-sea  interface  (<  few  wavelengths),  both  the  characteristic  flatness 
in  grazing  angle  and  any  obvious  resonance  behavior  within  a  given  frequency  band  can  disappear  (Figs.  12(a) 
and  12(c)).  This  suggests  an  alternate  explanation  for  the  enhanced  levels  (relative  to  rough  surface-scattering 
theory  predictions)  of  low-frequency  surface-scattering  strengths  that  have  been  observed  in  low-to-moderate 
sea  states.  Ecstatically,  the  scattering  picture  is  similar  for  fish  away  from  the  surface  (Fig.  12(f)),  and  more 
complicated  for  fish  very  near  the  surface  (Fig.  12(e)),  demonstrating  the  need  for  measuring  fish  behavior. 

In  general,  when  fish  are  near  the  air-sea  interface,  a  rich  variety  of  frequency  and  grazing-angle  behaviors 
of  <5™^  can  be  observed.  Figure  13  illustrates  this  with  a  real-world  example  of  backscattering  from  near¬ 
surface  salmon  during  a  low-sea-state,  night-time  measurement  in  the  Gulf  of  Alaska  during  the  Active  Adjunct 
to  Undersea  Surveillance  (AAUS)  Leg  5  Area  1  experiment  (May  1984)  [45].  These  displays  of  gSjqTAL 
demonstrate  that,  in  some  cases,  fish  can  be  the  dominant  scattering  mechanism  at  low  grazing  angles,  even  at 
low  densities— average  salmon  densities  in  this  environment  are  typically  a  few  hundred  individuals  per  square 
kilometer  [42].  (In  this  case,  with  by  up  to  20  dB  at  10°. )  It  further  illustrates  the  utility  of  the 

physics-based  scattering  model  as  an  inversion  tool.  After  inversion  for  the  salmon’s  density  and  depths  (in  this 
case,  2.5  to  7  m),  the  model  allows  extrapolation  in  frequency  and  angle.  In  this  case,  the  prediction  is  low 
backscatter  levels  at  very  low  grazing  angles  (<~5°). 

Fish  Near  the  Ocean  Bottom 

It  has  just  been  demonstrated  that,  in  some  cases,  backscatter  from  fish  near  the  sea  surface  can  dominate 
(and  be  confused  with)  surface  backscatter  from  the  air- water  interface  (and  bubble  clouds).  Similarly,  fish  near 
the  ocean  bottom  are  potentially  a  competing  scattering  mechanism,  both  to  the  bottom  interface  and  to  the 
sediment  volume.  Analogous  to  the  ocean-surface  case,  the  total  scattering  strength  will  be  modeled  as  the  sum 
of  contributions  from  the  ocean-bottom  interface  (o-^f ),  fish  (Oy^),  and  the  ocean-bottom  volume  (a^J): 

SS&TAL  =  10  •  log10  (ogr  +  +  a&) .  (38) 


(Here,  we  take  obv°0\  =  0 .) 

Figure  14  shows  a  low-frequency  example  of  backscattering  due  to  shoaling  rockfish  over  a  mudstone 
bottom.  This  instructive  example  derives  from  the  modeling  of  low-grazing-angle  bottom  backscattering  data 
collected  by  Kunz  of  NRL  during  LWAD  99-3  [28].  In  this  case,  the  acoustic  data  followed  a  psin0  angular 
dependence  (Fig.  14(a)),  an  empirical  relationship  commonly  ascribed  to  bottom  backscattering.  However, 
upon  modeling  the  scattering  contributions  of  both  the  rockfish  and  the  bottom,  it  was  found  [28]  that  this 
dependence  was  really  the  coincidental  combination  of  rockfish  scattering  at  low  angles  and  mudstone  scatter¬ 
ing  at  higher  angles. 
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Fig.  12  —  Monostatic  (a)-(d)  and  bistatic  (e)-(f)  predictions  of  the  angular  dependence  of  the 
scattering  strength  of  near-surface  salmon  (at  typical  Gulf  of  Alaska  sizes  and  densities  [42]) 
as  a  function  of  frequency  (150  to  1050  Hz)  for  two  sets  of  fish- layer  depths.  Subplots  (c)-(d) 
respectively  include  the  contribution  of  interface  scattering  (for  a  wind  speed  of  5  m/s)  to  the 
near-surface  salmon  scattering  displayed  in  (a)-(b). 
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Fig.  13  —  Data-model  comparison  for  near-surface 
salmon  +  interface  scattering  in  the  Gulf  of  Alaska. 
The  wind  speed  was  3.0  m/s.  The  acoustic  data 
were  collected  quasi-monostatically  using  omnidi¬ 
rectional  explosive  charges  as  sources,  which  deto¬ 
nated  at  depths  of  -457  m  below  a  horizontal  re¬ 
ceiver  towed  at  a  depth  of  -210  m.  For  reference, 
the  top  plot  includes  the  predicted  contribution  of 
interface  scattering  at  110  Hz  (dashed  curve). 


FREQUENCY  (Hz) 


Fig.  14  —  Predicted  (a) -(b)  angular  and  (c)  fre¬ 
quency  dependence  of  the  scattering  strength  of  rock- 
fish  uniformly  distributed  0-15  m  above  a  mudstone 
bottom  in  91  m  of  water.  Included  in  (a)-(c)  as  lightly 
dashed  curves,  are  the  free-field  rockfish  scattering 
responses  (i.e.,  not  including  bottom-interference 
effects).  Included  in  (a)  are  data  (symbols)  from 
LWAD  99-3,  [28,46]  and  a  thin  solid  curve  repre¬ 
senting  -27.5  +  10  •  loglosin0  dB.  In  (a)-(b),  the 
p’s  indicate  the  compressional-speed  critical  angles. 
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To  illustrate  this  sort  of  interplay  more  clearly,  Fig.  14  presents  a  simplified  version  of  the  LWAD  99-3 
scattering  scenario.  (Refs.  28  and  55  provide  more  details.)  Here,  we  consider  three  varieties  of  rockfish  (small, 
medium,  and  large)  in  91  m  of  water  that  are  uniformly  distributed  0  to  15  m  above  the  mudstone  bottom 

parameterized  in  Table  2.  (Here,  c0  =  1485  m/s.)  The  average  rockfish  layer  density  was  assumed  to  be  0.06 

fish/m2  (i.e.,  p ylsh  =  0.06  fish/m2  -s-  15  m  =  0.004  fish/m3).  Figure  14(a)  displays  monostatic  predictions  of  the 
backscattering  strength  at  500  Hz  from  the  rockfish  (thick  dashed  curve),  the  mudstone  bottom  (dotted  curve), 
and  their  total  (thick  solid  curve).  A  psinO  curve  matched  to  the  acoustic  data  (symbols)  overlays  these  curves. 
In  this  case,  we  see  that  which  mechanism  dominates  is  a  strong  function  of  grazing  angle  and  that,  even  at 
these  relatively  low  densities,  the  rockfish  dominate  the  scattering  below  15°.  At  increasingly  higher  densities, 
rockfish  scattering  would  dominate  the  scattering  over  an  increasingly  broader  range  of  angles. 

Ecstatically,  we  see  a  similar  picture.  Figure  14(b)  shows  that  the  basically  flat  dependence  of  rockfish 
scattering  strength  with  scattering  angle  implies  their  contribution  at  500  Hz  will  be  primarily  at  low  scattering 
angles  (and  minimal  near  specular).  Also  apparent  in  Figs.  14(a)  and  14(b)  are  local  peaks  in  the  rockfish’s 
scattering  response  associated  with  the  compressional-speed  critical  angle  (e.g.,  at  -49°). 

Another  factor  affecting  the  relative  dominance  of  fish  and  bottom  scattering  mechanisms  is  frequency. 
Figures  14(a)  and  14(b)  presented  predictions  for  500  Hz  that  correspond  to  the  peak  acoustic  response  of  the 
rockfish  at  these  depths.  Off  resonance,  a  fish’s  target  strength  is  reduced  (Fig.  14(c)),  so  that  which  scattering 
mechanism  dominates  can  depend  on  a  complex  mix  of  factors.  This  example  reinforces  the  usefulness  of 

physics-based  scattering  models  for  both  isolating  the  contributions  of  different  scattering  mechanisms,  and 

allowing  extrapolation  with  confidence  in  frequency  and  angle  (especially  bistatically). 

As  demonstrated  above,  the  acoustic  response  of  fish  can  be  mistaken  for  bottom  scattering  in  some  cases 
(especially  if  only  narrow  ranges  of  frequencies  or  angles  are  sampled).  In  particular,  since  fish  generate 
characteristic  broadband  frequency  responses,  their  scatter  could  be  confused  with  sediment-volume  scatter, 
also  generally  a  broadband  frequency-dependent  phenomenon  (a  cautionary  tale  for  those  doing  geo-acoustic 
inversions).  Furthermore,  for  a  nonuniform  distribution  of  fish,  say  shoaling  fish,  the  resultant  spatial  variabil¬ 
ity  in  could  be  mistakenly  attributed  to  sediment-volume  heterogeneities,  and  local  “hot  spots”  of  fish  to 
“geo-clutter.”  Figure  15  provides  a  data  example  from  October  2000  of  highly  localized  fish  schools  on  the 
New  Jersey  continental  shelf.  These  sonar  images  were  generated  by  Turgut  of  NRL  [56]  using  a  downward¬ 
looking,  2  to  8  kHz  chirp  echo  sounder,  which  pinged  twice  a  second  as  it  passed  over  the  bottom  at  a  rate  of 
~3  m/s.  (In  this  case,  the  backscatter  levels  from  the  fish  were  within  10  dB  of  those  from  the  sandy  bottom.) 


Fig.  15  —  Daytime  chirp  sonar  (2-8  kHz) 
record  displaying  fish  schools  adjacent 
to  the  ocean  bottom  on  the  New  Jersey 
Shelf  (~39°20'N,  72°50'W).  The  nomi¬ 
nal  water  depth  was  80  m  (courtesy  of  A. 
Turgut,  NRL  [56]) 
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Comments 

For  both  boundaries,  one  consequence  of  these  boundary-interference  models  is  that  they  predict  that  a 
fish’s  target  strength  is  generally  enhanced.  Typically,  the  effect  is  to  increase  the  average  scattered  intensity  by 
a  factor  of  4  (6  dB).  (When  Qinc  exactly  equals  Qscat,  this  factor  is  6  (~7.8  dB);  however,  as  these  two  angles 
begin  to  differ,  the  factor  quickly  converges  to  a  value  of  4.  Figure  14  provides  examples  of  this  enhancement— 
the  differences  of  the  thick  and  thin  dashed  curves.) 

Fish  near  the  ocean  surface  or  bottom  in  shallow  water  often  are  in  schools.  Modeling  and  experiments 
[4,57]  suggest  that  resonant  coupling  and  multiple  scattering  from  the  schools  can  cause  shifts  in  the  amplitude 
and  frequency  of  the  swimbladder  resonance.  This  could  result,  for  example,  in  increases  in  scattering  at 
frequencies  somewhat  below  an  individual  fish’s  expected  resonance.  So,  when  fish  school  in  the  presence  of  a 
boundary,  further  modification  to  the  above  formulas  will  be  required  to  accurately  model  their  scattering 
response. 

Because  fish  have  characteristic  frequency  dependencies  to  their  acoustic  responses,  broadband  measure¬ 
ments  combined  with  these  models  can  help  discern  the  extent  to  which  the  observed  scattering  is  due  to  fish 
or  to  the  ocean  boundary.  Also,  day/night  measurements  illustrating  fish  behavior  can  help  a  great  deal  in  this 
regard. 

6.  SURFACE  SCATTERING 

The  scattering  picture  in  the  dynamic  air-sea  boundary  zone  is  complex.  When  wave  breaking  is  signifi¬ 
cant,  scattering  from  the  rough  air-sea  interface  is  augmented  with  scattering  from  subsurface  bubbles.  At  high 
frequencies,  resonant  scattering  from  individual  bubbles  is  known  to  be  an  important  scattering  mechanism 
[58].  At  lower  frequencies  (<~5  kHz),  however,  observations  on  the  temporal  response  of  scattering  can  be 
used  to  rule  out  resonant  scattering  from  being  a  significant  scattering  mechanism.  In  this  case,  acoustic 
scattering  from  bubbles  can  still  occur  through  scattering  from  assemblages  of  bubbles  within  bubble  clouds,  or 
from  entire  clouds.  This  is  primarily  due  to  Bragg  scattering  from  inhomogeneities  in  the  sound-speed  field, 
which  is  proportional  to  the  Fourier  transform  of  the  index  of  refraction  field.  At  these  frequencies,  the  index 
of  refraction  depends  on  the  air-void  fraction,  but  not  on  the  details  of  the  bubble  distribution.  Hence,  a  key 
environmental  descriptor  of  bubble  clouds  at  these  lower  frequencies  is  the  air-void  fraction. 

Recent  low-frequency  (<  1500  Hz)  surface-scattering  measurements  under  the  CST  program  have  demon¬ 
strated  with  high  confidence  that  low-frequency  scattering  from  bubble  clouds  of  low  air- void  fraction  ( <  0.01  %) 
(microbubble  clouds)  is  responsible  for  the  high  reverberation  levels  observed  in  high  sea  states  [45,59].  In 
conjunction  with  the  key  acoustic  measurements,  Farmer  and  Vagle  of  the  Institute  of  Ocean  Sciences  made 
measurements  of  the  in  situ  air-void  fraction  that  they  have  been  able  to  relate  empirically  to  the  wind  speed 
[54] .  Spurred  by  these  CST  scattering  and  oceanographic  measurements,  NRL  has  developed  a  semi-empirical 
model  for  estimating  the  scattering  strength  due  to  bubble  clouds  that  is  parameterized  solely  by  the  wind  speed. 

This  bubble-cloud  scattering  model  has  been  integrated  with  our  model  for  into  a  broadband,  fully 
bistatic  surface  scattering  strength  (SSS)  model: 


SSS  =  10  •  log10  +  abub).  (39) 

The  resulting  model  is  parameterized  by  frequency  (10  to  5000  Hz),  the  scattering  angles  (incident,  scattered, 
bistatic),  and  environmental  variables  (wind  speed  and  two  surface-wave  spectral  parameters).  (As  will  be 
discussed  shortly,  a  further  approximation  allows  SSS  to  be  parameterized  solely  by  wind  speed.)  A  particular 
advantage  of  this  physics-based  model  is  its  reliance  on  readily  observable  oceanographic  variables  [60-61], 
quantities  that  can  be  measured  in  near-real-time  on  regional  scales  using  airborne  [62]  and  spaceborne  [63] 
remote  sensing. 
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This  model  offers  significant  advances  over  previous  (purely  empirical,  monostatic)  algorithms  [64-65], 
namely  physics-based  descriptions  of  scattering  from  both  the  rough  air-sea  interface  and  subsurface  bubbles; 
fully  bistatic  calculations;  and  applicability  to  a  broader  range  of  frequencies. 

In  this  section,  after  a  brief  summary  of  data  sources  used  in  deriving  the  semi-empirical  fits,  the  previ¬ 
ously  presented  model  is  illustrated  and  then  further  refined.  This  is  followed  by  a  description  of  our 
bistatic  model  for  abub.  These  models  are  then  used  to  briefly  explore  the  sensitivity  of  surface  scattering 
strength  to  environmental  variables  as  functions  of  the  acoustic  frequency  and  scattering  angles. 

Data  Sources 

The  deep-water  data  used  in  the  nonlinear,  multiparameter  curve  fitting  of  the  SSS  model  derive  from  a 
variety  of  sources:  Bachmann;  Brown  et  al. ;  Chapman  and  Harris;  Chapman  and  Scott;  Crowther;  and  Ogden 
and  Erskine— see  the  references  of  Refs.  53  and  58.  (The  curve-fitting  is  done  using  dB  values.) 

For  low  frequencies  (<1500  Hz),  the  primary  data  were  broadband  SUS  data  from  CST-4  and  CST-7 
Phase  2,  two  experiments  conducted  in  the  Gulf  of  Alaska.  In  both  of  them,  in  situ  measurements  of  bubble- 
cloud  properties  were  made,  and  the  fish  contributions  to  the  scattering  are  considered  to  be  well  understood 
[45]. 

Interface  Scattering  Model 

A  bistatic  scattering  model  for  was  described  above.  Here,  we  illustrate  its  sensitivities  to  its  inputs, 
and  derive  a  simplified  model  that  depends  environmentally  only  on  the  wind  speed  (and  c0). 

Model  Illustrations 

Figure  16  presents  representative  frequency  and  grazing-angle  behaviors  of  that  hold  for  most  wind 
speeds. 

Figures  16(a)  and  16(b)  illustrate  the  generally  moderate  frequency  dependence  of  air-sea  interface  scatter¬ 
ing.  In  each  plot,  the  three  sets  of  curves  correspond  to  grazing  angles  of  10°,  30°,  and  60°  for  a  wind  speed  of  20 
m/s.  (Similar  behavior  was  seen  at  most  wind  speeds.)  Figure  16(a)  shows  that,  as  expected  from  the  definition 
of  W(Qh ),  the  strongest  dependence  is  for  y2  =  3.0,  with  levels  rising  ~20  dB  as  the  frequency  increases  from 
100  to  10000  Hz,  and  the  weakest  dependence  is  for  y2  =  3.9  (~1  dB  increase).  This  is  reinforced  in  Fig.  16(b), 
where  y2  is  fixed  at  3.9  and  As  is  allowed  to  vary.  (Recall  the  only  frequency  dependence  in  our  model 
comes  from  y2.)  Figure  16(b)  also  shows  a  basically  linear,  monotonically  increasing  dependence  on  As  at  low 
angles,  with  a  more  complex  behavior  at  high  angles. 

Figures  16(c)  and  16(d)  show  the  very  strong  dependence  of  on  grazing  angle,  with  levels  increasing 
by  50  dB  or  more  as  the  grazing  angle  is  increased  from  5°  to  60°.  (This  behavior  holds  over  a  wide  range  of 
environmental  values  and  frequencies.) 

Re  the  grazing-angle  sensitivity  of  to  y2  and  As ,  Figs.  16(c)  and  16(d)  show  that  increasing  the  value 
of  each  parameter  leads  to  increasing  at  lower  grazing  angles,  while  at  higher  angles,  the  picture  is  more 
complicated.  Near  specular,  increasing  the  value  of  each  parameter  in  fact  leads  to  decreasing  .  (For 
scattering  from  rough  surface,  these  angular  trends  are  not  unexpected.)  Figure  16(a)  further  shows  that  for  a 
given  As,  the  angular  regions  where  there  will  be  a  monotonically  increasing  or  decreasing  dependence  of 
on  y2  will  be  frequency-dependent. 

Simplified  Model 

With  the  goal  of  constructing  a  SSS  model  that  depends  environmentally  on  only  the  wind  speed  U  (and 
c0),  this  model  was  fit  to  low  sea-state  CST-4  and  CST-7  Phase  2  data  to  derive  estimates  of  the  model’s 
two  surface-wave  spectral  parameters.  The  resulting  best- fit  values  were  y2  =  3.8  and  As  =  1.9  x  10-4  m3-s. 
(Initial  comparisons  with  in  situ  wave  spectral  values  have  been  favorable  [60-61].) 
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Fig.  16  —  Monostatic  predictions  of  interface  backscattering  strength  at  U  =  20  m/s  as 
a  function  of  (a)-(b)  frequency  and  (c)-(d)  grazing  angle,  for  two  sets  of  surface  wave 
spectral  variables.  In  (a)  and  (c),  As  =  0.0002  m3-s,  and  in  (b)  and  (d),  y2  =  3.9. 


Given  that  the  y2  derived  from  fitting  the  acoustic  data  was  3.8,  Fig.  16  suggests  that  will  generally 
have  a  mild  frequency  dependence  (at  least  in  the  open  ocean).  This  is  illustrated  in  Fig.  17,  which  provides  a 
representative  example  from  CST-7  Phase  2  of  broadband,  low- wind-speed  scattering  strengths  observed  during 
the  CST  series  of  experiments. 

Bubble-Cloud  Scattering  Model 

Figure  18  provides  a  low-frequency  example  of  how  the  surface  backscattering  strength  can  change  when 
wave  breaking  is  significant.  In  this  high-wind-speed  example  from  CST  7  Phase  2,  we  see  the  strongly 
enhanced  scattering  levels  relative  to  those  predicted  as  being  due  to  the  air-sea  interface  (black  curves).  We 
also  note  that  in  this  case,  the  maximum  response  is  not  at  the  highest  frequency,  but  at  962  Hz,  and  that,  even 
at  82  Hz,  the  presence  of  bubbles  contributes  significantly  to  the  backscattering.  We  now  describe  a  bistatic  abub 
model  that  attempts  to  capture  these  features. 

The  bubble-cloud  scattering  prediction  derives  from  the  stochastic  model  of  Gilbert  [66].  In  Ref.  66,  the 
mean  bubble-cloud  backscattering  cross  section  is  modeled  as  a  product  of  an  effective  horizontal  wavenumber 
spectrum  and  a  “geometric  factor.”  (However,  as  Gilbert  notes,  very  little  is  known  about  the  former.)  This 
geometric  factor  describes  how  the  backscatter  intensity  from  a  bubble  cloud  depends  on  two  competing  factors: 

•  The  rapid  decrease  of  average  air-void  fraction  with  depth,  which  suggests  that  most  scattering  occurs  very 
close  to  the  surface;  and 
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Fig.  17  —  Low- wind- speed  surface  scattering  strengths 
measured  in  the  Gulf  of  Alaska  during  CST-7  Phase  2 
(February  1992).  The  wind  speed  was  4.5  m/s.  The 
acoustic  data  were  collected  quasi-monostatically  us¬ 
ing  omnidirectional  explosive  charges  as  sources,  which 
detonated  at  depths  of  -540  m  below  a  horizontal  re¬ 
ceiver  towed  at  a  depth  of  -220  m.  Overlaying  the  data 
are  bistatic  SSS  model  predictions  (withT5  =  0.00019 
m3-s,  y2  =  3.8)  at  83  Hz  (thick  black,  solid  curve)  and 
1350  Hz  (thick  black,  dashed  curve). 


Fig.  18  —  High- wind- speed  surface  scattering  strengths 
measured  in  the  Gulf  of  Alaska  during  CST-7  Phase  2 
(February  1992).  The  wind  speed  was  17.9  m/s.  The 
acoustic  data  were  collected  quasi-monostatically  using 
omnidirectional  explosive  charges  as  sources,  which  deto¬ 
nated  at  depths  of  -560  m  below  a  horizontal  receiver 
towed  at  a  depth  of  -260  m.  Also  plotted  are  bistatic 
model  predictions  (with  As  =  0.00019  m3-s,  y2  = 
3.8)  at  82  Hz  (thick  black,  solid  curve)  and  1378  Hz 
(thick  black,  dashed  curve). 


•  The  vanishing  of  the  total  acoustic  field  at  the  surface  and  its  increase  to  a  maximum  several  meters  below 
the  surface,  because  of  interference  between  the  incident  and  surface-reflected  waves  (Lloyd’s  mirror 
pattern). 

Assuming  that  a  bubble  cloud  may  be  modeled  as  a  vertical  distribution  of  uncorrelated  point  scatterers  (at 
depths  z),  the  convolution  of  these  two  effects  in  the  monostatic  case  leads  to  a  total  intensity  of  [66] 

/  =  \  sin4 ( k0z  sin 6) exp(-2z  / d)dz  =  A2k*d5 /[(l  +  k^d2)(l  +  4k^d2)]  (40) 

where  d  is  the  e-folding  depth  of  the  air  fraction,  kv  =  k0  sin  0  is  the  magnitude  of  the  vertical  wavenumbers  of 
both  the  incident  and  scattered  fields,  and  the  At  are  scaling  factors. 

The  approach  we  have  taken  is  to  use  the  above  geometric  factor  as  a  starting  point  for  a  new  bistatic  ahub 

algorithm.  First,  Eq.  (40)  was  extended  to  3D  geometries  [67].  Then,  the  scaling  factor  was  taken  to  be  a  three- 

parameter  function  of  k0  and  d.  To  determine  these  three  ohub  parameters,  we  needed  a  ready  means  of 
estimating  d.  Fortunately,  there  is  a  prescription  relating  d  to  the  wind  speed  U.  Based  on  measurements  during 
CST-7  Phase  2,  Farmer  and  Vagle  have  derived  the  following  empirical  formula  [54]: 

d  =  0.557  -  0.117  *  17  +  0.0109  •  U2  U  >  7.5  m/s, 

=  -  0.19509  +  0.06503  •  U  3  <  U  <  7.5  m/s,  (41) 

=  0  U  <  3  m/s. 

Using  this  formula,  multidimensional,  nonlinear  regression  routines  were  then  used  to  fit  our  SSS  model  to  the 
full  data  set  (over  10000  data  points).  (Because  no  in  situ  measurements  of  d  were  made  during  the  non-CST- 
7  Phase-2  measurements,  Eq.  (41)  was  assumed  to  hold  for  all  the  acoustic  data.)  The  result  is  [59] 
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_  0.0019J51  V-6^2,^[6  +  3(^,  +  +  (k^  -kj,)2d4] 

bub  2(1  +  k2vid2 )(1  +  k2sd2)[l  +  (kvi  -kVj)2d2][l  +  (kvi  +  kv/)2d2\ 

where  kv  i  =  -k0  sin  Qinc  and  kvs  =  -k0  sin  Qscat.  As  with  the  3D  and  Gh/-xh  models,  the  point-scatterer  and 
flat-interface  assumptions  imply  that,  to  first  order,  abub  depends  on  Qinc  and  §scaV  but  not  §bi. 

Thus,  for  this  model,  abub  depends  on  two  environmental  parameters:  d  and  c0,  which,  via  a  further 
approximation  (Eq.  (41)),  becomes  just  U  and  c0. 

Comments 

For  practical  applications  of  this  model  as  kernels  in  reverberation  or  performance  models,  a  field  of 
discrete  bubble  clouds  is  treated  as  a  spatially  (and  temporally)  uniform  layer,  which  in  turn  is  treated  as  an 
effective  modifier  of  surface  properties.  (Being  a  semi-empirical  model,  the  nonuniform  nature  of  real  bubble 
fields  is  already  embedded  in  the  derived  values  of  the  abub  parameters.)  Hence,  it  can  be  used  as  a  surface 
scattering  strength  with  no  need  to  introduce  a  separate  layer  into  the  modeled  waveguide. 

Model-Model  Comparisons 

The  residual  standard  error  rse  for  fitting  the  resultant  new  SSS  model  to  the  full  deep-water  data  set  was 
2.4  dB,  compared  to  4.1  dB  when  applying  the  Chapman-Harris  [64]  empirical  formula.  (When  restricted  to 
data  corresponding  to  frequencies  below  1000  Hz,  the  new  model  had  an  rse  of  2.2  dB,  compared  to  4. 1  and  3.1 
dB  when  applying  the  Chapman-Harris  and  Ogden-Nicholas-Erskine  [65]  empirical  formulas,  respectively.) 

Simplified  SSS  Model 

We  have  just  derived  a  SSS  model  that  relies  environmentally  on  just  U  and  cQ .  However,  as  the  dependence 
on  c0  g  [1400,  1600]  m/s  is  typically  weak,  it  could  be  hardwired  at  1500  m/s,  with  only  minor  (<~1  dB) 
effects  on  SSS.  Doing  so,  the  result  is  a  simplified  SSS  model  that  depends  environmentally  on  only  the  wind 
speed. 

Model  Illustrations 

With  Fig.  19  we  begin  to  look  at  the  relative  contributions  of  bubble  clouds  and  the  interface  to  surface 
scattering  strength.  Figures  19  through  22  present  monostatic  predictions,  while  Fig.  23  presents  bistatic 
predictions.  For  the  predictions  of  interface  scattering,  we  set  As  =  0.0002  m3-s  and  y2  =  3.9  since,  of  the 
parameter  values  of  Fig.  16,  they  are  the  most  representative  of  open-ocean  values.  (In  Figs.  19  through  23, 
there  are  no  bubble-cloud  scattering  predictions  for  wind  speeds  less  than  3  m/s,  as  our  model  assumes  obub  = 
0  for  U  <  3  m/s.  For  all  predictions,  we  take  c0  =  1500  m/s  and  §bi  =  180°.) 

Grazing-Angle  Dependence 

Figure  19  displays  the  grazing-angle  dependence  at  1500  Hz  for  five  wind  speeds  (2.5,  5,  10,  15,  and  20 
m/s)  of  bubble-cloud  scattering  strength  (Fig.  19(a)),  interface  scattering  strength  (Fig.  19(b)),  and  the  sum  of 
the  two,  SSS  (Fig.  19(c)).  Figure  19(a)  illustrates  the  relatively  flat  grazing  angle  dependence  of  bubble-cloud 
scattering,  and  that  the  range  of  grazing  angles  over  which  the  flatness  occurs  increases  with  increasing  U.  It 
also  shows  the  especially  strong  dependence  of  abub  on  U  at  low  grazing  angles  (<~30°  at  1500  Hz). 

In  contrast,  g™?  has  both  a  milder  and  more  complex  dependence  on  U  as  a  function  of  grazing  angle 
(Fig.  19(b)),  with  increasing  with  U  at  the  lower  angles  and  decreasing  with  U  near  specular  (as  would  be 
expected  as  the  surface  becomes  increasingly  rougher). 

The  combination  of  effects  (Fig.  19(c))  reinforces  the  general  rule  of  thumb  that  bubble  clouds  dominate 
the  scattering  at  low  grazing  angles  (when  wave  breaking  is  significant)  and  the  rough  air-sea  interface  at  high 
angles,  with  a  wind- speed-dependent  mix  at  moderate  angles. 
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GRAZING  ANGLE  (deg) 


Fig.  19  —  Monostatic  predictions  of  surface 
backscattering  strength  at  1500  Hz  as  a  func¬ 
tion  of  grazing  angle  for  five  wind  speeds,  (a) 
Bubble-cloud  backscattering  strength;  (b)  and 
(c)  correspond  to  interface  and  total  (bubble- 
cloud  +  interface)  backscattering  strength,  re¬ 
spectively,  assuming  As  =  0.0002  m3-s  and  y2 
=  3.9. 


FREQUENCY  (Hz) 

Fig.  20  —  Monostatic  predictions  of  surface 
backscattering  strength  at  0  =  10°  as  a  func¬ 
tion  of  frequency  for  five  wind  speeds,  (a) 
Bubble-cloud  backscattering  strength;  (b)  and 
(c)  correspond  to  interface  and  total  (bubble- 
cloud  +  interface)  backscattering  strength,  re¬ 
spectively,  assuming  As  =  0.0002  m3-s  and  y2 
=  3.9. 
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WIND  SPEED  (m/s) 

Fig.  21  —  Monostatic  predictions  of  surface 
backscattering  strength  at  1500  Hz  as  a  func¬ 
tion  of  wind  speed  for  five  grazing  angles,  (a) 
Bubble-cloud  backscattering  strength;  (b)  and 
(c)  correspond  to  interface  and  total  (bubble- 
cloud  +  interface)  backscattering  strength,  re¬ 
spectively,  assuming  As  =  0.0002  m3-s  and  y2  = 
3.9. 


0  5  10  15  20 

WIND  SPEED  (m/s) 

Fig.  22  —  Monostatic  predictions  of  surface 
backscattering  strength  at  three  grazing  angles 
as  a  function  of  wind  speed  for  five  frequen¬ 
cies.  For  comparison,  monostatic  predictions 
of  just  the  contributions  due  to  the  air-sea  in¬ 
terface  at  50  and  5000  Hz  (thin  solid  curves) 
are  included  in  each  subplot.  (The  interface¬ 
scattering  predictions  assumed  As  =  0.0002 
m3-s  and  y2  =  3.9.) 
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Fig.  23  —  Modeled  bistatic  surface  scattering  strengths  at  1500  Hz  as  a  function  of 
scattered  grazing  angle  Qscat  for  two  incident  grazing  angles.  In  (a)-(b),  bubble-cloud 
scattering  strength  is  shown  (solid  curves)  for  four  wind  speeds  (U  =  5,  10,  15,  and 
20  m/s)  and  interface  scattering  strength  (dashed  curves)  for  two  winds  (U  =  5  and  20 
m/s),  assuming  As  =  0.0002  m3-s  and  y2  =  3.9.  In  (c)-(d),  the  total  (bubble-cloud  + 
interface)  scattering  strength  is  shown  for  five  wind  speeds.  (The  dashed  vertical  lines 
indicate  the  specular  directions  in  each  case.) 


Frequency  Dependence 

Figure  20  shows  the  dependence  of  SSS  at  low  grazing  angles  on  frequency  and  wind  speed.  As  in  the 
previous  figure,  Fig.  20(a)  corresponds  to  bubble-cloud  scattering,  Fig.  20(b)  to  interface  scattering,  and  Fig. 
20(c)  to  the  sum  of  the  two  (SSS).  Again,  the  set  of  five  curves  in  each  plot  corresponds  to  five  wind  speeds 
(2.5,  5,  10,  15,  and  20  m/s).  A  grazing  angle  of  10°  is  picked  as  an  interesting  one  where  both  bubble  clouds 
and  the  interface  can  have  a  significant  influence  on  the  acoustics. 

Figure  20(a)  illustrates  the  strong  frequency-dependence  of  bubble-cloud  scattering  below  1  kHz,  with  a 
peak  occurring  at  a  frequency  that  decreases  with  increasing  wind  speed.  (Both  points  are  dramatically  illus¬ 
trated  in  Fig.  18,  with  peak  acoustic  response  in  this  case  at  962  Hz.)  Above  1  kHz,  the  frequency  dependence 
is  relatively  mild  for  wind  speeds  greater  than  10  m/s.  Figure  20(b)  shows  that,  for  our  given  wave  spectrum, 
the  frequency  and  wind  speed  dependence  of  interface  scattering  are  fairly  benign.  (For  a  different  wave 
spectrum,  the  wind- speed  dependence  would  be  similar  to  those  in  Fig.  20(b),  although  the  curves  would 
exhibit  a  different  dependence  on  frequency  (Fig.  16(a))  and/or  be  at  displaced  levels  (Fig.  16(b).) 

The  effects  of  including  both  scattering  mechanisms  are  shown  in  Fig.  20(c).  At  this  grazing  angle,  at  wind 
speeds  above  10  m/s,  bubble-cloud  scattering  dominates  above  -250  Hz.  Below  10  m/s,  surface  scattering 
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strength  can  be  due  to  a  combination  of  bubble-cloud  and  interface  scattering,  while  below  5  m/s,  interface 
scattering  dominates  (except  at  high  frequencies).  As  the  grazing  angle  increases,  interface  scattering  becomes 
increasingly  dominant  (Fig.  19). 

Wind-Speed  Dependence 

Figure  21  illustrates  the  wind-speed  dependence  of  SSS  as  a  function  of  grazing  angle.  As  in  the  previous 
figure,  Fig.  21(a)  corresponds  to  bubble-cloud  scattering,  Fig.  21(b)  to  interface  scattering,  and  Fig.  21(c)  to  the 
sum  of  the  two  (SSS).  In  this  case,  the  set  of  five  curves  in  each  plot  corresponds  to  five  grazing  angles  (5,  10, 
15,  30,  and  60  deg).  A  frequency  of  1500  Hz  is  picked  as  an  interesting  one  in  which  both  bubble  clouds  and 
the  interface  can  have  a  significant  influence  on  the  acoustics. 

Figure  21(a)  illustrates  the  very  strong  wind- speed  dependence  of  bubble-cloud  scattering  strength,  espe¬ 
cially  in  the  regime  from  when  wave  breaking  becomes  significant  (U  ~  5  m/s  for  1500  Hz)  until  the  onset  of 
high  sea  states  (U  ~  15  m/s  for  1500  Hz),  where  abub  begins  to  appreciably  slow  its  increase.  In  Fig.  21(b), 
equivalent  curves  are  presented  for  interface  scattering  strength;  these  reflect  our  previous  comments  as  to 
interface  scattering’s  strong  grazing-angle  and  relatively  mild  wind-speed  dependence  (compared  to  those  of 
bubble-cloud  scattering). 

Figure  21(c)  shows  the  general  dominance  of  interface  scattering  at  high  grazing  angles  and  that  of  bubble- 
cloud  scattering  at  low  grazing  angles  in  moderate-to-high  winds.  The  contrast  among  these  three  figures  shows 
the  complexity  of  surface  scattering  at  low-to-moderate  winds  and  low  grazing  angles;  the  frequency  depen¬ 
dence  of  obub  adds  to  the  fun. 

Figure  22  again  plots  the  dependence  of  SSS  on  wind  speed,  but  this  time  as  a  function  of  frequency  at  three 
low  grazing  angles.  This  figure  reinforces  the  points  that  bubble  clouds  become  an  increasingly  important  driver 
of  backscattering  strength  with  both  increasing  frequency  and  wind  speed,  and  with  decreasing  grazing  angle.  For 
example,  at  a  grazing  angle  of  5°,  the  contribution  of  bubble-cloud  backscattering  to  SSS  for  frequencies  >  ~1  kHz 
can  be  30  to  40  dB  higher  than  the  contribution  of  air-sea-interface  backscattering  for  wind  speeds  in  excess  of  15 
m/s. 

Bistatic-Angle  Dependence 

Figure  23  illustrates  that  the  general  trends  discussed  so  far  extend  to  bistatic  geometries.  In  these  examples 
at  1500  Hz,  bubble-cloud  and  interface  scattering  strengths  (Figs.  23(a-b)),  and  their  total  (SSS)  (Figs.  23(c-d)), 
are  shown  for  a  set  of  wind  speeds  for  two  incident  grazing  angles  of  10°  and  60°  (thus,  specular  is  at  170°  and 
120°,  respectively). 

A  comparison  of  Figs.  23(a)  and  23(c)  with  Fig.  19  shows  that,  in  general,  at  high  scattering  angles  (toward 
specular),  interface  scattering  dominates,  while  at  low  scattering  angles,  bubble-cloud  scattering  dominates  for 
wind  speeds  of  10  m/s  and  greater.  This  is  emphasized  in  Figs.  23(b)  and  23(d),  where  the  incident  grazing 
angle  has  been  increased  to  60°,  with  the  expected  result  that  interface  scattering  dominates  over  a  larger  range 
of  scattering  angles  and  wind  speeds. 

Comments 

We  have  presented  a  new  semi-empirical  broadband  surface  scattering  strength  formula  parameterized 
environmentally  by  the  wind  speed,  two  surface-wave  spectral  parameters,  and  the  air-void  fraction  e-folding 
depth  d.  Further  approximations  led  to  a  formula  whose  sole  environmental  parameter  is  the  wind  speed  U. 
From  an  analysis  of  the  variance  of  CST-4  and  CST-7  Phase  2  SSS’s  measured  during  periods  of  similar  wind 
speeds  (~9  to  11  m/s),  it  can  be  shown  that  wind  speed  is  not  the  whole  story.  (For  example,  the  relative  wind 
and  wave  directions  can  potentially  be  of  significance:  Farmer  and  Vagle  have  found  [54]  enhanced  wave 
breaking  (bubble  generation)  when  the  wind  direction  changes,  and  reduced  wave  breaking  when  the  wind  and 
wave  directions  are  opposed.)  Hence,  to  provide  more  accurate  predictions  of  abub ,  in  situ  measurements  of  d 
and  other  environmental  measures  of  the  air-sea  boundary  conditions  will  be  necessary.  Furthermore,  it  should 
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be  noted  that  the  relationship  between  d  and  U  can  depend  on  the  site  and  the  season  [54] ,  so  that  physics-based 
versions  of  this  relationship  need  to  be  developed. 

The  SSS  model  presented  should  be  viewed  as  an  interim  model.  For  example,  we  are  considering  different 
roughness  spectra  for  the  model,  and  are  trying  other  functional  forms  for  sbub.  We  also  need  to  properly 
include  the  coherent  component  in  the  model  to  predict  forward  scattering  (surface  loss).  (Our  new  model 
for  dispersed-fish  scattering  in  the  presence  of  the  sea  surface  is  contributing  to  the  ongoing  effort  to  definitively 
identify  the  relative  contributions  of  bubbles,  waves,  and  fish  in  low-to-moderate  sea  state  data  from  the  CST 
and  other  archival  data  sets.) 

Subtle  differences  in  low-  and  mid-frequency  SSS  characteristics  (such  as  in  grazing  angle  behavior)  were 
observed  when  separately  deriving  fits  for  data  above  and  below  1500  Hz.  These  phenomena  (perhaps  sugges¬ 
tive  of  different  scattering  physics)  are  under  study. 

7.  CONCLUSIONS 

These  new  bistatic  scattering  strength  models  should  help  advance  the  accuracy  and  robustness  of  predic¬ 
tions  of  multistatic  active  sonar  performance.  By  having  a  physics  basis,  the  models  allow  extrapolation  in 
frequency  and  to  most  3D  geometries.  Hence,  they  should  be  particularly  useful  for  shallow-water  applications 
where  a  wide  variety  of  scattering  angles  is  encountered,  but  where  making  controlled  measurements  of 
scattering  strength  is  difficult.  Our  models’  broadband  and  bistatic  natures  permit  their  application  to  a  wide 
range  of  multistatic  active  sonar  systems  and  concepts. 

The  models  have  proved  essential  for  isolating  scattering  mechanisms  and  thereby  advancing  our  under¬ 
standing  of  the  complex  acoustic  interaction  processes  at  the  ocean  boundaries.  (For  example,  being  able  to 
systematically  identify  and  remove  fish-scatter  contributions  from  surface-scatter  data  has,  in  turn,  led  to  more 
accurate  surface-scattering  models.) 

Because  of  their  basis  in  physical  theory,  these  models  can  be  used  as  inversion  (“remote-sensing”)  tools. 
(For  example,  fits  of  the  interface  surface  scattering  model  to  broadband,  low-wind-speed  acoustic  data  were 
used  to  back  out  representative  surface-wave  spectral  values.)  For  algorithm  development  and  in  situ  operations, 
these  models  have  several  potential  uses: 

•  Benchmarking  geo-acoustic  inversion  and  environmentally  adaptive  (EA)  algorithms 

•  Forming  the  basis  for  “next-generation”  EA  techniques 

•  Extrapolation  in  frequency  and  angle  (especially  bistatically),  given  an  environmental  parameter  estimation 
from  a  small  band  of  angles  at  a  few  frequencies 

•  Helping  to  invert  for  hard-to-estimate  geophysical  parameters  (such  as  the  spectral  strength  and  exponent 
of  bottom  roughness) 

•  Identifying  what  oceanographic  and  geophysical  quantities  need  to  be  measured,  and  to  what  accuracy 

•  Serving  as  a  guide  to  expected  acoustic  behavior  (e.g.,  identifying  when  elastic  effects  are  significant) 

•  Serving  as  a  guide  in  asset  choice  and  deployment,  and  frequency  and  waveform  selection  (optimizing 
signal-to-interference  ratios). 

Despite  their  apparent  complexity  in  some  cases,  the  models  are  straightforward  to  code  and  run  quite  quickly 
(especially  if  vectorized— we  have  MATLAB®  and/or  S-PLUS®  codes  for  all  models). 

In  general,  measures  of  sonar  performance  depend  nonlinearly  on  the  reverberation,  which  in  turn  depends 
nonlinearly  on  the  environmental  (oceanographic,  geophysical,  and  biological)  variables.  A  key  benefit  of  these 
models  is  that  they  allow  a  systematic  determination  of  the  relative  influence  of  these  environmental  inputs  on 
the  strength  of  the  acoustic  scattering.  In  turn,  used  as  scattering  submodels  in  active  sonar  performance 
models,  they  allow  a  more  accurate  estimation  of  the  relative  influence  of  the  environment  on  sonar  perfor¬ 
mance. 

Furthermore,  by  independently  varying  the  values  of  the  environmental  parameters,  the  resultant  impact  on 
the  scattering  can  be  estimated  in  a  statistical  sense.  Hence,  the  relative  variability  or  uncertainty  in  sonar 
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performance  can  be  assessed,  and  the  expected  variance  modeled  (i.e.,  helping  to  estimate  reverberation 
probability  density  functions). 

We  now  provide  a  categorized  list  of  the  major  points  and  recommendations  of  the  paper: 

General  Interface  Scattering 

Bottom  and  surface  interface  scattering  were  well-modeled  by  the  SSA.  Its  use  is  strongly  recommended 
over  empirical  (such  as  Lambert’s  or  Lommel-Seeliger’s  law  for  the  bottom)  and  PT  approaches. 

To  more  accurately  model  scattering  strength  at  high  scattering  angles  and  in  the  specular  direction,  the 
pure  power-law  wave  spectrum  needs  to  be  modified  to  include  a  low-wavenumber  cut-off.  Furthermore,  for  the 
sea  surface,  more  general,  directional  water-wave  spectra  need  to  be  incorporated  into  the  scattering  model. 
(These  efforts  are  underway  at  NRL.) 

Bottom  Interface  Scattering 

Scattering  from  the  bottom  interface  has  a  complex  and  strong  dependence  on  the  incident  and  scattered 
grazing  angles,  and  on  the  geophysics.  Of  the  latter,  at  low  grazing  angles  the  spectral  parameters,  sounds 
speeds  and  density  will  typically  have  the  most  influence  on  the  acoustics.  (Generally,  the  harder  the  bottom, 
the  stronger  the  acoustic  response  at  low  scattering  angles.)  The  dependence  of  the  scattering  response  on 
acoustic  frequency  is  generally  monotonic,  with  levels  typically  increasing  moderately  with  increasing  fre¬ 
quency. 

Elastic  effects  can  be  significant  below  the  critical  angle  for  bottom-interface  scattering  from  shear¬ 
supporting  bottoms,  i.e.,  skin  effects  due  to  evanescent  shear  waves  can  be  significant.  Surface  roughness 
enhances  these  effects.  (Our  SSA  model  has  the  advantage  of  accurately  predicting  such  effects  in  cases  where 
the  acoustic  wavelength  cannot  resolve  violations  of  the  single-valuedness  assumption  implicit  in  the  “rough- 
surface”  assumption.)  One  particular  manifestation  is  shear-induced  nulls  in  scattering  strength  at  low  grazing 
angles.  As  illustrated  with  the  mudstone  example  from  the  Oregon  coast,  these  nulls  can  occur  even  when  the 
shear  speed  is  well  below  the  water  sound  speed. 

Volume  Scattering 

Bladdered  fish  can  be  significant  scatterers  at  low  grazing  angles,  even  at  relatively  low  densities.  When 
they  are  near  the  ocean  surface  or  bottom,  their  scattered  intensity  can  be  significantly  altered  (by  a  factor 
between  0  and  16;  typically  4).  This  response  depends  strongly  on  the  incident  and  scattered  grazing  angles,  the 
distance  of  the  scatterers  from  the  boundary,  the  acoustic  frequency,  and  for  the  ocean  bottom,  its  physical 
properties.  Due  to  the  more  reflective  nature  of  the  ocean  surface,  near-surface  fish  will  elicit  a  more  variable 
response  than  near-bottom  fish. 

Near-ocean-surface  fish  can  be  a  significant  and  complex  low-  and  mid-frequency  scattering  mechanism. 
For  example,  they,  and  not  bubbles,  may  in  fact  be  responsible  for  reported  deviation  of  low-frequency  surface 
scattering  strengths  at  low  grazing  angles  from  rough-surface  scattering  predictions  in  low  and  moderate  sea 
states. 

Similarly,  near-ocean-bottom  fish  could  be  responsible,  in  some  cases,  for  scattering  attributed  to  interface 
and/or  sediment- volume  scattering.  (Recall  the  above  scattering  strength  data  with  a  psinG  dependence  that 
suggested  a  coincidental  combination  of  rockfish  scattering  at  low  grazing  angles,  and  mudstone  bottom  scatter¬ 
ing  at  intermediate  angles.)  As  was  demonstrated,  bladdered  fish  very  near  the  bottom  (or  even  resting  on  it) 
can  still  be  appreciable  scatterers.  Furthermore,  because  of  their  often  high  spatial  variability  in  shallow  water, 
local  “hot  spots”  of  near-bottom  fish  could  be  mistaken  for  “geoclutter” 

In  general,  fish  will  be  of  more  importance  as  a  source  of  interference  to  active  sonar  systems  at  mid¬ 
frequencies  and  in  shallow  water.  The  strength  of  volume  scattering  will  be  strongly  dependent  on  local 
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conditions  (such  as  the  types,  number  densities,  sizes,  and  depths  of  fish  present,  and  the  time  of  day). 
Furthermore,  their  spatial  and  temporal  variability  will  in  turn  lead  to  acoustic  variability. 

Surface  Scattering 

A  surface  scattering  strength  model  was  presented  that  offers  significant  advances  over  previous  (purely 
empirical,  monostatic)  algorithms  in  several  ways:  (a)  the  scattering  descriptions  of  both  the  rough  air-sea 
interface  and  subsurface  bubbles  derive  from  fundamental  physics;  (b)  the  algorithms  that  implement  them 
operate  in  3D  mode;  and  (c)  the  resulting  models  apply  over  a  broader  frequency  range  (which  includes  ASW 
frequencies).  (Furthermore,  for  scenarios  where  the  models  all  apply,  the  new  model  more  accurately  fits  the 
observed  data.)  While  this  model  represents  a  considerable  improvement  in  predictive  capabilities,  it  should  be 
noted  that  it  is  still  semi-empirical,  with  a  significant  amount  of  its  key  oceanographic  and  low-frequency 
acoustic  data  deriving  from  one  open-ocean  location  at  basically  one  time  of  the  year. 

Scattering  due  to  bubble  clouds  can  be  dominant  at  low  grazing  angles  when  wave  breaking  is  significant. 
Bubble  clouds  become  an  increasingly  important  driver  of  backscattering  strength  with  both  increasing  fre¬ 
quency  and  wind  speed,  and  with  decreasing  grazing  angle.  For  quantitative  predictions,  a  good  estimate  of  the 
in  situ  wind  speed  is  needed.  (Wind-history  effects  are  relatively  unimportant.)  Bubble  cloud  scattering  typi¬ 
cally  has  a  relatively  flat  response  with  grazing  angle  (especially  for  wind  speeds  above  10  m/s).  The  range  of 
angles  where  the  response  is  flat  increases  with  both  wind  speed  and  frequency.  Bubble-cloud  scattering  has  a 
strong  dependence  on  frequency  below  ~1  kHz  (with  a  relatively  flat  response  above  for  wind  speeds  in  excess 
of  10  m/s). 

Scattering  due  to  the  rough  air-sea  interface  depends  most  strongly  on  the  scattering  angle  and  environ¬ 
mental  conditions.  To  predict  its  strength,  a  good  estimate  of  the  surface  wave  spectrum,  primarily  at  short 
wavelengths,  is  needed.  At  large  scattering  angles,  the  interface  is  the  dominant  surface-scattering  mechanism 
at  all  frequencies,  even  in  high  sea  states.  Interface  scattering  has  a  relatively  mild  dependence  on  environmen¬ 
tal  parameters.  (ohub  will  be  much  more  responsive  to  locally  changing  oceanographic  conditions,  such  as  an 
increase  in  wind  speed,  than  will  .) 

A  key  feature  of  the  surface  scattering  strength  model  presented  is  its  environmental  reliance  on  readily 
observable  oceanographic  variables,  which  in  turn  offers  the  potential  for  accurate  predictions  of  SSS  in  near- 
real-time  on  regional  scales  through  the  use  of  remote  sensing.  (The  accuracy  and  availability  of  real-time, 
wide-area  measures  of  these  variables  are  continually  improving  [62-63].) 

For  most  frequencies  and  sea  conditions,  the  model  results  demonstrated  that  at  low  scattering  angles,  the 
key  environmental  quantity  to  be  measured  is  the  wind  speed.  At  moderate-to-high  scattering  angles,  bubble- 
cloud  scattering  is  less  important.  In  this  regime,  both  the  wind  speed  and  directional  wave  spectra  need  to  be 
measured. 

Issues 

While  these  models  provide  improved  predictions  of  mean  (incoherent)  scattering  levels,  technical  issues 
remain.  A  primary  need  is  high-quality  acoustic  and  environmental  data  to  provide  the  ground  truth  necessary 
to  rigorously  evaluate  the  models,  and  to  assess  the  generality  and  limitations  of  their  physical  assumptions.  For 
each  interface,  the  coherent  (specular)  component  needs  to  be  properly  combined  with  the  incoherent  compo¬ 
nent  to  accurately  predict  boundary  loss.  Specific  issues  include: 

Bottom  Scattering 

We  have  explored  in  this  paper  how  the  properties  of  the  bottom  drive  the  acoustics.  To  improve  the 
accuracy  and  robustness  of  bottom- scattering  models,  much  work  remains: 

•  Wide  ranges  of  measured  acoustic  and  geophysical  values  have  been  observed,  even  in  environments  of 

similar  “bottom  type”  (e.g.,  sand).  Bottom  properties  are  often  difficult  to  measure  in  situ  to  the  level  of 
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detail  required  by  models— needed  are  new  measurement  techniques  to  characterize  the  bottom.  This,  in 
particular,  includes  measuring  interface  roughness. 

•  Consequently,  scattering-model  validation  at  low-  and  mid-frequencies  is  poor,  principally  due  to  the  lack 
of  high-quality,  concurrent  acoustic  and  geophysical  measurements.  (There  is  a  particular  sparsity  of  3D 
acoustic  data.) 

•  The  statistics  of  sediment-volume  scattering  is  not  well  understood— needed  are  3D  geologic  models  of 
sediment  heterogeneities. 

•  How  significant  is  subcritical  angle  penetration  due  to  surface  roughness  [68]  at  low-  and  mid-frequencies? 

•  For  acoustic  measurements,  it  is  difficult  in  shallow  water  to  isolate  direct-path  scattering  (particularly  at 
low  frequencies  and  angles)  [36].  Range  dependence  makes  it  even  tougher. 

Volume  Scattering 

To  predict  the  scattering  contributions  from  fish  in  a  given  environment,  knowledge  of  fish  behavior  as  well 
as  fish  types  and  densities  are  required.  (For  example,  knowing  the  diurnal  depth  behaviors  and  seasonal 
migration  patterns  of  particular  fish  species  can  aid  greatly  both  in  predicting  their  acoustic  behavior  on 
regional  scales,  and  in  mitigating  their  influence  on  sonar  performance.)  Particular  issues  include: 

•  Modeling  the  dynamism  of  fish  and  fish  schools  in  shallow  water  to  whatever  accuracy  is  realistically 
possible. 

•  Developing  the  capability  for  at-sea  assessments  of  fish  school  statistics. 

•  Developing  and  validating  fish  school  scattering  models  that  incorporate  effects  due  to  the  presence  of  the 
ocean  boundaries. 

Surface  Scattering 

To  provide  more  accurate  predictions  of  scattering  from  bubble  clouds,  in  situ  measurements  of  the  air- void 
fraction  and  other  environmental  measures  of  the  air-sea  boundary  conditions  will  be  necessary.  Hence,  the  air- 
void  fraction  represents  a  candidate  quantity  for  estimation  with  remote-sensing  technologies  (perhaps  using 
lidar).  The  relationship(s)  between  the  air- void  fraction  and  the  wind  speed  (and  other  environmental  descrip¬ 
tors)  can  be  site  and  seasonally  dependent  [54] ,  so  that  physics-based  models  of  the  Bragg-scale  structure  of  the 
void  fraction  need  to  be  developed. 

To  provide  more  accurate  measurements  of  scattering  from  the  air-sea  interface,  in  situ  estimates  of  2D 
surface- wave  power  spectra  are  needed  (e.g.,  via  synthetic  aperture  radar  (SAR)  imagery). 

A  deeper  understanding  of  bubble  properties  near  the  surface  is  needed.  Important  outstanding  questions 
include: 

•  What  properties  of  bubbles  drive  the  acoustics? 

•  How  are  they  affected  by  environmental  conditions  such  as  temperature  and  salinity? 

•  How  do  they  evolve  temporally?  Spatially? 

•  At  what  frequencies  are  individual  resonant  bubbles  important?  Clouds  of  bubbles?  What  are  the  relative 
roles  of  attenuation  and  Bragg  scattering  generated  by  heterogeneities  in  the  sound-speed  profile  as  a 
function  of  frequency  and  sea  state? 

•  What  are  the  differences  between  the  open  ocean,  shallow  water,  and  the  surf  zone?  (In  coastal  waters,  SSS 
levels  are  often  significantly  enhanced  [58].) 

•  How  do  surface  waves  couple  into  the  problem? 

•  What  is  the  relationship  between  scattering  and  noise? 
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Recommendations 

We  close  with  some  recommendations  for  any  scattering  measurement: 

•  Maximize  the  frequency  and  grazing-angle  (both  incident  and  scattered,  if  possible)  coverage  to  sort  out 
scattering  mechanisms  and  help  invert  for  environmental  parameters. 

•  Perform  day /night  measurements  to  help  sort  out  the  fish  contributions. 
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